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Abstract. In this paper we formulate necessary conditions for the integrability in the Jacobi 
sense of Newton equations q = —F(q), where q € C™ and all components of F are polynomial and 
homogeneous of the same degree /. These conditions are derived from an analysis of the differential 
Galois group of the variational equations along special particular solutions of the Newton equations. 
We show that, taking all admissible particular solutions, we restrict considerably the set of Newton's 
equations satisfying the necessary conditions for the integrability. Moreover, we apply the obtained 
conditions for a detailed analysis of the Newton equations with two degrees of freedom (i.e., n = 2). 
We demonstrate the strength of the obtained results analyzing general cases with deg Fi = I < 4. 
For I = 3 we found an integrable case when the Newton equations have two polynomial first 
integrals and both of them arc of degree four in the momenta p\ — qi, and P2 — Q2- Moreover, 
for an arbitrary k, we found a family of Newton equations depending on one parameter A. For an 
arbitrary value of A one quadratic in the momenta first integral exist. We distinguish infinitely 
many values of A for which the system is integrable or superintegrable with additional polynomial 
first integrals which seemingly can be of an arbitrarily high degree with respect to the momenta. 

1 Introduction 

In this paper we study the integrability of the following class of Newton's equations 

q = -F(q), q = ( qi ,...,q n ) T G C\ (1.1) 

In spite of the fact that such equations already appeared at the beginning of the mechanics 
in the Newton Principia [28j |. there is not too many rigorous results concerning their 
integrability. The only exception is the case of potential forces, i.e., the case when 

dV 

FM = ^ to), 



for a certain function V. In fact, for potential forces equations (II. ip are equivalent to 
Hamilton's equations 

dV 

Q = P, P = -^(q), (1-2) 

with natural Hamiltonian 

H = ^p T p + V(q). (1.3) 

First of all, for Hamiltonian systems we have a precise notion of the integrability i.e. the 
integrability in the Liouville sense. Moreover, for such systems we have powerfull methods 
and tools for integrability studies. This is why in this paper we investigate the integrability 
problem for Newton's equations with non-potential forces. However, for such equations 
'integrability' is not well defined, so we should first specify what, in the considered context, 
the integrability means. 

The oldest notion of the integrability is related to quadratures. We say that the 
considered system of differential equations is integrable by quadratures if all its solutions 
can be obtained in a closed form by a finite number of algebraic operations (including 
inversion of functions) and quadratures, i.e. integrals of 'known' functions, see e.g. [3]. 
In fact, already in his Principia Newton reduced to quadratures the simplest problems of 
dynamics (see theorems XXXIX, XLI, LIII and so on, see also comments in 0]). A similar 
definition of integrability was applied by Poincare [29l ]. 

The integrability by quadratures was a prototype for other definitions of integrability. 
As a matter of fact, many of these new definitions can be considered as necessary conditions 
for the integrability by quadratures. Generally, if for a given dynamical system 

x = v(x), x = (xi, . . .x n ) T £ C n , (1.4) 

there exists a sufficient number of tensor invariants, i.e. tensor fields that are constant 
along phase curves, then the system is solvable by quadratures. The most known two 
examples are the following. If system (ll.4p possesses n — 1 functionally independent first 
integrals, then it is integrable by quadratures. The Lie theorem says that if system (|1 .4[) 
has n symmeties, i.e., if there exist vector fields u\{x) = v(x),U2(x), . . . , u n {x) which are 
linearly independent and [ui,Uj] = 0, for i,j = 1, . . . ,re, then system (jl.4p is integrable 
by quadratures. 

In [lj| C. G. J. Jacobi developed, introduced by L. Euler, a method of the integrable 
factor. The Jacobi Last Multiplier of equation (|1.4p is an invariant ra-form 

fM = h(x)dxi A • • • A dx n . 

Invariance of \x means that 

n d 

dw(hv) : =^_(^) = 0. 

i=l * 

Theorem 1.1 (Jacobi). Assume that system (II. 4p possesses n—2 functionally independent 
first integrals and a Last Jacobi Multiplier. Then it is integrable by quadratures. 
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We can also recall the well known fact that the integrability of a Hamiltonian system 
in the Liouville sense implies its integrability by quadratures, see e.g. [lj. For more details 
about the integrability by quadratures see [2|;lZ|]. 

For the purpose of this paper we introduce the following definition. 

Definition 1. We say that system f 1 1 . 4 [) is integrable in the Jacobi sense if it admits an 
invariant n-form and possesses n — 2 functionally independent first integrals. 

This definition of the integrability is very often applied for nonholonomic systems, see 

e.g. EBHSQ. 

Let us notice that Newton's equations (II. ip can be written in the following form 

q=p, p=-F(q), (1.5) 

so they admit 

H = dqi A • • • A dq n A dpi A • • • A dp n , 

as an invariant 2n-form. Thus, in a case when equations (jl.5|) are not Hamiltonian, we 
say that they are integrable if they are integrable in the Jacobi sense. 

The integrability does not depend on the chosen coordinates, however, the form of 
equations (|1.5j) and the homogeneity of F do. A linear change of variables given by 
q = AQ, p = AP, where A G GL(n, C), transforms equations (jl.5p into 

Q = P, P = —Fa(Q), where F A (Q) := A~ 1 F(AQ). (1.6) 

Hence it is reasonable to divide all forces into disjoint equivalence classes. We say that 
forces F and G are equivalent iff there exists A G GL(ra, C) such that G = Fa- Later, 
specifying a force we just give a representative of equivalent forces in the above sense. 

Our main purpose is to consider Newton's equations which are not Hamiltonian. Hav- 
ing this in mind, we introduce the following definition. 

Definition 2. We say that a force F is potential iff there exists A G GL(n, C) and 
V G C[q] such that F A = W. 

For Hamiltonian systems there are many approaches to the integrability problem: the 
Hamilton-Jacobi theory, perturbation techniques, normal forms, splitting of separatrices 
and, recently, the Morales-Ramis theory based on differential Galois theory. Many in- 
tegrable Hamiltonian systems were found and exhaustively analyzed by means of these 
methods. For non-Hamiltonian systems we have not so many various approaches. 

The main aim of this paper is to formulate necessary conditions of integrability in the 
Jacobi sense for Newton systems ()1 . 5[) . In this aim, we use an approach which is based 
on an investigation of variational equations for a particular complexified solution. This 
approach was developed by by S.L. Ziglin (351 : I3H ] for Hamiltonian systems. However, 
the basic idea of his method can be used for a study of the integrability of an arbitrary 
dynamical system. The fact that the considered system possesses a meromorphic first 
integral implies that there exists a rational function which is an invariant of the monodromy 
group of the variational equations. This gives a restriction on the monodromy group. 
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For the first time this method was applied for proving the non-integrability of a non- 
Hamiltonian system by S.L. Ziglin in 37; 3^; 40; 39]. 



The Ziglin theory appeared to be a very strong tool for proving the non-integrability 
of Hamiltonian systems. Thanks to works of A. Baider, R. .C. Churchill, J. J. Morales, 
J. -P. Ramis, D. L. Rod, C. Simo and M. F. Singer the Ziglin theory was considerably 



developed, see [J; |24|; [25|; |26|] and references therein. The main idea of this development 
is following. As in the Ziglin theory, we investigate the variational equations along a non- 
equilibrium solution. But, instead of the monodromy group of the variational equations, 
we investigate their differential Galois group. This change gives the substantial profit 
because it is much easier to investigate the differential Galois group of the given equations 
than their monodromy group, and the differential Galois group contains the monodromy 
group as a subgroup. The main theorem due to Morales and Ramis states that if a 
Hamiltonian systems is integrable in the Liouville sense, then the identity component of 
the differential Galois group of variational equations along a non-equilibrium solution is 
Abelian. Some parts of the described differential Galois approach to the integrability can 
be adopted to non-Hamiltonian systems. The key implication is the following. If a system 
possesses a meromorphic first integral, then the differential Galois group of variational 
equations along a non-equilibrium solution has a rational invariant. In [19J the reader 
finds the first application of the differential Galois approach for proving non-integrability 
of a non-Hamiltonian system. This type of the integrability analysis was also applied 
for a nonholonomic system called the Suslov system in [21]. In this paper we apply the 
differential Galois approach for integrability studies of Newton's equations. 

To formulate the main results of this paper we have to fix our assumptions. Hence, 
if it is not otherwise stated, we assume that force F = (iq, . . . ,F n ) is polynomial and 
homogeneous, i.e., Fi G C[q], and degFj = / > 1 for i = 1, . . . , n. Moreover, we always 
write deg-Fj = I = k — 1, and letter k is used only in this context (if F is a potential force, 
then k is the degree of the potential). Thus we have always k > 2. By a direction in C n 
we understand an equivalence class of parallel non-zero vectors. Of course a direction is 
determined just by one non-zero vector. 

Definition 3. A direction d £ C n is a Darboux point of a force F(q), iff d is parallel to 
F(d), i.e., d A F{d) = 0, and F{d) + 0. 

To find all Darboux points of a given force it is enough to find all non-zero solutions 
of equations 

F(d) = d, (1.7) 

and take those which give different directions. To simplify exposition, we always use such 
normalization that a Darboux point d of F satisfies the above equation. The set of all 
Darboux points of a given force F is denoted by T>f- It can be an empty, infinite or finite 
set. If T>f is finite, then for a generic F it has D(n, k) := [(k — l) n — l]/(k — 2) elements. 
A Darboux point d defines a particular solution of system (|1.5p by 



q(t) = <p(t)d, p(t) = ip(t)d, (li 
where ip(t) satisfies the equation 



<p = -ip k \ 
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<p 2 = y(e-<p k ), (1.9) 



This equation determines a family of hyperelliptic curves 

2 

k 

depending on a parameter e £ C*. 

The variational equations along solution (II. 8p have the form 

x = y, y = -<p(t) k - 2 F'(d)x, (1.10) 

or simply 

x = - v (t) k - 2 F'{d)x, (1.11) 

where F'(d) is the Jacobi matrix of F calculated at a Darboux point d. Let us assume 
that this matrix is diagonalizable. Then, in an appropriate basis equations (|1.11|) have 
the form 

Vi = ->*<p(t) k ~ 2 m, i = l,...,n, (1.12) 

where Ai, . . . , A n are eigenvalues of F'{d). By homogeneity of F one of the eigenvalues, 
let us say A ra , is (k — 2). We denote the rest of eigenvalues by A = X(d) = (Ai, . . . , A n _i) G 

Our first theorem is the following. 

Theorem 1.2. Assume that the Newton system (jl.ip with polynomial homogeneous right- 
hand sides of degree greater than one is integrable in the Jacobi sense. If d is a Darboux 
point of F, such that F'(d) is semi-simple, then the identity component of the differential 
Galois group of variational equations along the particular solution defined by d is Abelian. 

Our next theorem gives a computable criterion for the integrability in the Jacobi sense. 

Theorem 1.3. Assume that the Newton system (|1.1|) with polynomial homogeneous right- 
hand sides of degree I = k — 1 > 1 is integrable in the Jacobi sense. If d is a Darboux 
point of F, such that F'{d) is semi-simple, and Ai, . . . , A n are eigenvalues of F'(d), then 
(k, Ai) for i = 1, . . . , n belong to the following list 

case k A 



1. 


k 


k 

P + -jPiP ~ 1) 






2. 


k 








3. 


3 




1 3 
~24 + 32 


(1 + 4p) 2 








1 3 
~24 + 50 


(2 + hp) 2 


4. 


4 








5. 


5 




9 t 1 
~40 + 10 


(2 + hp) 2 



(1.13) 
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where p is an integer. Moreover, among Ai, . . . A n at most two belong to the first item of 
the above table. 

The above theorem follows from Theorem 1 1 . 2 1 and the fact that each of the second order 
variational equations in (|1.12|) can be transformed into a Gauss hypergeometric equation 
for which the differential Galois group is known. 

Our next result is of a different nature. It is easy to understand that an application 
of Theorem 11.31 to a general multi-parameter family of F gives rather an unsatisfactory 
result: we distinguish infinitely many families of F that satisfy conditions of Theorem [L3j 
One can try to find more than one Darboux point and apply Theorem 11.31 to all of them. 
But in practice, without any guiding idea, an analysis of the obtained conditions is at least 
difficult. 

To describe our approach, let us define for a Darboux point d of F quantities Aj = 
Aj — 1, for i = l,...,n — 1, where Aj are eigenvalues of F'(d). We denote them by 
A = A(d) = (Ai, . . . , A n _i) G C™ -1 . Let t%, for i = 0, ...,n — 1, denote elementary 
symmetric polynomials in (n — 1) variables. The following theorem shows that for a given 
F, quantities A(d) calculated at different Darboux points are not arbitrary. It gives n 
relations among A(d) calculated over all Darboux points d £ Dp. 

Theorem 1.4. Assume that F has exactly D(n,k) Darboux points d £ T>p. Then A(d) 
satisfy the following relations: 

E ^r^^^r^+k-iy, (i.i4) 

or, alternatively 

T r(A(d)) = fn-i- 1 

r n _ x (A(d)) 1 ] ^ 

d£V F 8=0 

for < r < n — 1. 



The most important consequence of this theorem is the following. We show that for a 
fixed k, the set of solutions of (|1.14p or (|1.15p . such that for each d S T>f ah components 
of \(d) satisfy conditions of Theorem 11.31 is finite. Moreover, there is an algorithmic way 
to find all of those solutions. In other words, the above theorem gives a possibility to 
investigate the integrability of a general polynomial Newton system with a fixed degree of 
homogeneity. 

It is worth mentioning here that the main idea which allows to derive the result con- 
tained in Theorem 11.41 is to consider the following auxiliary system 

= F(q) (1.16) 
associated with Newton's equations (]1.5p . It is a homogeneous first order system, so we 



can perform the Kovalevskaya analysis for it, see e.g. 17]. Then it appears that the 
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introduced quantities A, are the Kovalevskaya exponents calculated at Darboux points of 
system ()1.16j) . For more details, see Section [5j 

Without doubt the integrability is not a generic phenomenon. Thus one would like 
to weaken the assumptions of Theorem 11.41 and consider forces F with a smaller number 
of Darboux points. We do not know how to treat such cases for an arbitrary n, but for 
n = 2 we present an approach which allows to formulate a theorem much stronger than 
Theorem 11.41 for details see Section 

Applications of our general theorems for Newton's equations with two degrees of free- 
dom enabled us to distinguish for an arbitrary k a one parameter family for which the 
necessary integrability conditions formulated in Theorem 1 1 . 3 1 are seemingly sufficient. An 
analysis given in Section [7.41 supports the following conjecture. 

Conjecture 1.1. Assume that (k,X) belongs to an item of table f|1.13j) . Then Newton 
equations 

Qi = Pi, Pi = -A<7i<?2~ 2 92 =P2, P2 = q^ 1 (1-17) 
are integrable in the Jacobi sense with polynomial first integrals I\ and I2, where 

1 2 1 k 
h = - 2 Pl + pl 

Moreover, for an arbitrary M > we find A such that the degree of I2 with respect to 
momenta is greater than M and there is no an additional polynomial first integral inde- 
pendent with I\ and degree with respect momenta smaller or equal M . Additionally, if 
(k,\) belongs to an item different from 1 in table (|1.13p . then there exist two additional 
polynomial first integrals I2 and I3 which are functionally independent together with I\ . 

The plan of this paper is the following. In Section [2] definitions and some general results 
concerning the relations between first integrals of dynamical systems and invariants of their 
differential Galois group are presented. In Section [3] particular solutions and variational 
equations are calculated. Section H] is crucial because contains necessary conditions of 
the Jacobi integrability. In Section [5] consequences of the existence of more particular 
solutions related to Darboux points for the integrability analysis are considered. Section [6] 
is a compendium of all general results for systems of two Newton's equations. In the 
next section we apply these results to the general classes of two Newton's homogeneous 
equations with the right-hand sides of degree two and three and for such systems the 
integrability analysis is almost complete. Some partial general results for two Newton 
equations with higher degrees of homogeneity are also presented. Appendix A contains 
two approaches to the problem how to recognized that the analyzed Newton's system is 
Hamiltonian. In Appendix B the basic facts about higher order variational equations and 
their applications are shown. The necessary integrability conditions expressed by means 
of higher order variational equations appeared to be a very effective tool for the analysis of 
Hamilton equations. In this part of the paper we formulate a conjecture that they can be 
also applied for the Newton equations and examples show the validity of this hypothesis. 
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2 Basic facts from general theory 



Let us consider a system of differential equations 

—x = v(x), x£UcC n , teC, (2.1) 
at 

where v is a holomorphic vector field in the considered domain U. For a non-equilibrium 
particular solution <p(t) we consider also the variational equations along (pit), i.e., 

^ = A(t)*i, A(t) = ^{t)). (2.2) 

The considered particular solution <p(t) defines a Riemann surface T immersed in C n with 
t as a local parameter. The right hand sides of variational equations (|2.2p depend in fact 
on a point V. 

If F : C n D W — » C is a holomorphic function defined in a certain connected neigh- 
borhood W of the solution <p(t), then, by definition, the leading term / of F is the lowest 
order non-vanishing term of expansion 

n<p® + o = + o(ii£ir +1 )> F m ^ o, (2.3) 

i.e., /(£) := .F m (£). Note that /(£) is a homogeneous polynomial of variables £ = 
(£i,...,£ n ) of degree m; its coefficients are polynomials in (pit). If F is a meromor- 
phic function, then F = P/Q for certain holomorphic functions P and Q. In this case, the 
leading term / of F is defined as / = p/q, where p and q are leading terms of P and Q, 
respectively. In this case /(£) is a homogeneous rational function of To express this fact 
more precisely we denote by IK = M(r) the field of meromorphic functions on V. Then / 
is a rational function of £i, . . . , £ n with coefficients in K, i.e., / G K(£) := K(£i, . . . , £„). 

It is not difficult to prove that if F is a meromorphic first integral of equation (12.11) . 
then its leading term / is a first integral of variational equations (|2,2|) . Moreover, if 
Fi,...,Fk are functionally independent meromorphic first integrals of (|2.1|) . then, by 
the Ziglin Lemma, we can assume that their leading terms fu ■ ■ ■ , A are functionally 
independent first integrals of (|2.2p . For proofs and details see [33; 0; 3; 13]. 

Let S C GL(n,C) denote the differential Galois group of 1^. If f(g(£)) = f(£) for 
every g £ S, then we say that / is an invariant of S. The following lemma says that first 
integrals of variational equations are invariants of their differential Galois group. It was 
proved in [24[], see also Lemma III. 1.13 on p. 63 in [3J]. 

Lemma 2.1. If equation (|2.ip has k functionally independent first integrals which are 
meromorphic in a connected neighborhood of a non- equilibrium solution (p(t), then the 
differential Galois group 9 of the variational equations along pit) has k functionally in- 
dependent rational invariants. 

For a fixed i £ C (or, for a fixed point (pit) on the Riemann surface T), a first integral 
/ of the variational equations is an element of C(£). The set of all rational invariants of 
S is a field denoted by C(£) s . 



8 



The differential Galois group 9 of a system of linear equations is a linear algebraic 
group, see e.g. So, in particular, it is also a Lie group. This fact allows to 

reformulate the necessary conditions for the integrability in the language of Lie algebras. 
However, passing from a Lie group to its Lie algebra, we have to pay something. Namely, 
we take into account only the identity component 9° of 9- In particular, if 9 is finite, then 
the Lie algebra of 9 is trivial, i.e., it consists of one zero vector. 

Let q C GL(n, C) denote the Lie algebra of 9- An element Y £ q can be considered as 
a linear vector field: x i— ► Y(x) := Yx, for x £ C n . We say that / £ C(oj) is an integral 
of g, iff Y(/)(a;) = df(x) • y(ic) = 0, for all Y £ g. All rational first integrals of the Lie 
algebra q C GL(ra, C) form a field denoted by C(x) s . 

Proposition 2.1. If fi,...,fk £ C(a;) are algebraically independent invariants of an 
algebraic group 9 C GL(n, C), then they are algebraically independent first integrals of the 
Lie algebra q of 9- 

Proof. Let X £ q. Then g(t) = exp[tX] defines for t £ M, \t\ small enough, a one parameter 
subgroup of 9- Thus if / £ C(x) s , then 

f(exp[tX]x) = f(x), 

for any X £ g. Differentiating the above equality with respect to t at i = 0, we obtain 
X(f) = 0, i.e., / £ C(a!)« □ 

It can happen that two systems, one integrable and the other non-integrable, have 
the same variational equations along a chosen particular solution. In such situation, to 
show non-integrability of the second system we have to apply either another method or 
another particular solution or investigate higher order variational equations. The last 
approach appears very effective. In this paper we use it analyzing several examples and it 
is described in Appendix B. 

3 Differential Galois group of variational equations 

We showed in Introduction that the variational equations (|1 . 12j) for a particular solu- 
tion (|1.8p corresponding to a Darboux point d are a direct product of second order equa- 
tions of the form 

fj = -\<p(t) k -\ A£C, (3.1) 

where <p(t) satisfies (|1.9p . Let 9(&,A) C Sp(2, C) denote the differential Galois group of 
the above equation and S(k, A)° the identity component of $(k,\). Then the differential 
Galois group 9 of variational equations (|1.12p is a direct product 

9 = 5(k, Ai) x • ■ • x 9(k, A n ) C Sp(2n, C). (3.2) 

Hence, it is clear that in order to know the properties of 9, we should know as much as 
possible about S(k,X). Fortunately, the properties of the last group can be described in 
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details. Following Yoshida [34| it is convenient to define a new independent variable z 
in (J3ZED by 

z := ±<p(t) k , 

where (p{t) satisfies f|1.9f) . Then equation (|3.ip transforms into a linear second order 
equation with rational coefficients 

x „ (k - 1 3k - 2 \ . Xi , . 

This equation is a special case of the Gauss hypergeometric equation 

z(l - z)rj" + [c - (a + b + l)z]r]' - abrj = 0, (3.4) 

with 

a + 6=^, a6 = --, c=l--. (3.5) 
Thus, the differences of exponents at singularities z = 0,1 and oo for equation ()3.4p are 
1 ,1 .1 



p = 1 — c = - , a = c — a — b = -, r = a - b = — yf {k — 2) 2 + 8/cA, 

respectively. Let G(k,X) C GL(2,C) denote a differential Galois group of (13.31). The 
differential Galois group of the hypergeometric equation is well known, see e.g. 
However, to make this fact usefull, we should know a relation between 9(k, A) and G(k, A). 
It can be shown, see Proposition 4.7 in 0], that the identity components of these groups 
are the same. 

The well known Kimura theorem [lit ] specifies all cases when all solutions of the hy- 
pergeometric equation (|3.4p are Liouvillian. In these cases, by the Lie-Kolchin theorem, 
the identity component of the differential Galois group is solvable. A direct application of 
the Kimura theorem to equation (13. 3p gives the following. 



Lemma 3.1. The identity component G(k,X)° of the differential Galois group of equa- 
tion (|3.3p with k > 2 is solvable if and only if (k, X) belong to the list given in Theorem \1.3[ 

The above lemma is not sufficient for our purpose. We show the following. 

Lemma 3.2. Assume that (k,X) belongs to the list given in Theorem \1.3\ and k > 2. 
Then, for the differential Galois group G(k,X) of equation (|3.3p . the following statements 
hold. 

1. G{k, X)° is Abelian. 

2. If {k, X) belongs to the item 1 of the list, then 



G(k, A)° = Ti 



1 c 
1 



c€C . (3.6) 
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3. If (k, A) belongs to items 2-5 of the list, then G(k, A)° = {E}. 



Here E is the identity matrix. 
Proof. At first we transform equation (|3.3f) into the normal form. To this end we put 



w = rj exp 



- / jjdz 



._ c-(a + 6+l)z 



2 

where a, b and c are given by (|3.5p . Then we obtain 

„ _ P 2 -l + ^(l-P 2 -r 2 + ^ 2 ) + z 2 (r 2 -l) 
™ " 4z 2 (z-l) 2 ^ 

For this equation exponents at 0, 1 and at infinity are 

1(1-^1(1 + ^)1, 11(1-^1(1 + ^)1 |_l(i_ r) _l(i + r )l (3.9) 



respectively. Monodromy and differential Galois groups of (|3.8[) are now subgroups of 
SL(2,C). It is important to remark here that the identity components of the differential 
Galois group of (|3.3|) and (|3.8|) are the same. Notice also that the differences of exponents 
at singular points were unchanged. We denote by G(k, A) the differential Galois group of 
equation (|3.8|) . 

By assumption G(k, A)° is solvable, and thus G(k, A)° is also solvable. Suppose that 
G(k, A)° is solvable but not Abelian. Then, by Theorem 4.12 on p. 31 in 15J] and Propo- 
sition 4.2 in [32|, there is only one possibility: G(k,\) = G(k,X)° = T, where T is the 
triangular subgroup of SL(2, C). So, such case can appear only if the considered equation 
is reducible. Using the well known criterion for the reducibility of Riemann P equation, 
see e.g. 16j, we easily obtain that equation (|3.8p is reducible iff A = p + kp(p — 1) /2 where 

P ez. 

We show that if equation (|3.8p is reducible, then G(k, A) is a proper subgroup of 7, 
and in this way we prove the first statement of our lemma. 

First, let us notice that in the considered reducible case the respective exponents at 
singular points 0, 1 and infinity are following 

1 _ JL i . JL\ JIH [ 2 + fc(Z + 2) 2 + k(l-2) \ 

2 2fc'2 + 2fc/' \4'4J' \ 4k ' 4k J ' 1 ' 

where I is an odd integer. The difference of exponents at the singular point z = 1 is 1/2. 
Thus, by Lemma 4.3.6 on p. 90 in 12| equation (|3.8p has a solution of the form: 

w = z r (l- z) s h{z), 

where h{z) is a polynomial, and r is an exponent at z = 0, and s in an exponent at z = 1, 
As r and s arc rational, there exists j G N such that w J G C(z). Now, by Proposition 4.2 
in [2I], Q(k, A) is either a proper subgroup of the diagonal group, or a proper subgroup of 
the triangular group. This finishes proof of the first statement of the lemma. 
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To prove the second statement of we have to show that 9(k, A) cannot be a subgroup 
of the diagonal group. Assume that it is diagonal. Then, in particular, the monodromy 
group of equation f|3.8|> is diagonal. This last group is generated by two elements Mq, 
Mi G SL(2, C) — the monodromy matrices corresponding to canonical loops encircling 
singular points and 1, respectively. If the monodromy group of equation (13. 8p is diagonal, 



then by Lemma 4.3.5 on p. 90 in [12|], at least one of the matrices Mq, M\ or MqM± is 
a scalar matrix, i.e., ±E. Taking into account, that the eigenvalues of M, are exp2-7rirj ) i 
and exp27rirji, where rn and are exponents at z = i, we easily conclude that for 
k > 2 it is impossible. This finishes the proof of the second statement. 

To show the last statement it is enough to notice that if (k, A) belongs to item 2 of 
the list (|1.13p . then, using Theorem 2.9 (b) on p. 525 in [6J] one can check that G(k,\) 
is finite. Moreover, from the proof of the Kimura theorem given in [lf| it follows that 
if (k,X) belongs to items 3-5 of the list 11.131 then G(k, A), and thus G(k, A), is a finite 
primitive group. As the identity component of a finite group is just the identity element, 
this finishes the proof. □ 

Passing from group G(k, A) to its Lie algebra from the lemma proved above we have 
as an immediate consequence the following corollary. 

Corollary 3.1. Let g be the Lie algebra of the differential Galois group G(k,X) of equa- 
tion (|3.3|) . Assume that (k, A) belongs to the list given in Lemma \'3.1\ and k > 2. If (k, A) 
belongs to the item 1, then dining = 1, and otherwise dining = 0. 

4 Necessary conditions for the Jacobi integrability 

The main purpose of this section is to give proofs of Theorems 11.21 and 11.31 For this 
purpose we develop appropriate tools investigating certain Lie algebras. 

4.1 Certain Poisson algebra 

An element Y of Lie algebra sp(2n, C), considered as a linear vector field, is a Hamiltonian 
vector field given by a global Hamiltonian function H : <C 2n — ► C, which is a homoge- 
neous polynomial of 2ti variables • • • , x n , y±, . . . , y n ) of degree 2, i.e. H G C2[x,y] := 
C2[xi, . . . , x n , yi, . . . , y n \. In this way, we identify Lie algebra sp(2n, C) with a C-linear 
vector space C2[au, y] with the canonical Poisson bracket as the Lie bracket. Thus, for a 
Lie algebra q C sp(2n, C) ~ C2[x,y], a rational function / £ C[a;,y] is a first integral of 
g, iff {H, /} = 0, for all H G g. 

Now, we consider a more detailed case when q is a Lie subalgebra of sp(2,C). It is 
easy to show that Lie algebra sp(2, C) does not admit any non-constant first integral. 

Proposition 4.1. A rational function f G C(x,y) is a first integral o/sp(2,C), iff f G C. 

Proof. Let / G C(x, y) be a first integral of sp(2, C) ~ C2[x, y\. Thus, {/, H} = 0, for each 
H G C2[x,y]. Let us take H = x 2 . Then, 

{f > H] = - 2x % = °' 
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and this shows that / does not depend on y, i.e., / G C(x). Taking H = y 2 , we show that 
/ does not depend on x. Hence / G C. □ 

The above proposition shows that only proper subalgebras of sp(2, C) can have non- 
constant first integrals. 

Proposition 4.2. If g is a Lie subalgebra o/sp(2,C) and dimcfl > 0, then the number of 
algebraically independent rational first integrals of g is not greater than one. 

Proof. As dinicg > 0, there exists a non-zero H G sp(2,C) ~ C2[x, y\. The number of 
rational algebraically independent first integrals of a non-zero linear Hamiltonian vector 
field Xh in C 2 is at most one. □ 

Proposition 4.3. If g is a Lie subalgebra o/sp(2,C) and dimc0 = 2, then C(x,y) s = C. 

Proof. All two dimensional Lie algebras are solvable, so g is solvable. Thus a connected 
Lie group G C sp(2, C) with Lie algebra g is solvable. By the Lie-Kolchin theorem G is 
conjugate to the triangular group 

b 



7:-- 



a 

a' 1 



a G 



b G C 



The Lie algebra t of T is isomorphic to g, and is generated by two elements 

hi 



"1 " 


, h 2 = 


"0 


1" 


-1 









(4.1) 



(4.2) 



Let Hi and H 2 be Hamiltonian functions from C2[x,y] such that linear vector fields Xh 1 
and Xh 2 are represented by matrices hi and h 2 , respectively. It is easy to check that 



Hi = xy, 



H 2 = -y\ 



(4.3) 



We show that C(x, y) { = C. Assume that there exists / G C(x,y) t \C Hence {/, Hi} = 
for i = l,2. But 

{/,*>=»g=0, 
so, / G C(y). However, for / G C(y), we have 

{/ , ifl} = - 9 | =0 , 

and this implies that / G C. A contradiction with the assumption that / is not a constant 
shows that C(x,y) 1 = C. Moreover, as the Lie algebras t and g are isomorphic, we have 
also C{x,y) = C. □ 

Let K D C be a field. Then, for the Lie algebra sp(2, C) ~ C 2 [x, y], we can also consider 
first integrals which belong to K(x,y). Here we assume that 



da da 
— = and — 

ox oy 



for all a G K. 
For our further considerations we need the following lemma. 
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Lemma 4.1. Let g be a one dimensional Lie subalgebra 0/ sp(2, C) ~ C2[x, y\. If f G 
is a rational first integral of g, then there exists a nonzero element h G C[x,y] such 
that f G K(h). 

Proof. Let us assume that g is nilpotent, i.e., g is generated by H = y 2 . Then we have 

= {ff,/} = -2,g. 

Thus, / does not depend on x, so for this case we choose h = y. 

The only other possibility is that g is diagonal, i.e., it is generated by H = xy. First, 
let us assume that / is a polynomial in (x, y). We can write / uniquely as a sum of homo- 
geneous components. It is clear that if / is the first integral of H, then each homogeneous 
component of / is also the first integral of H. Thus, let us assume that / is homogeneous 
of degree s and let us represent it in the form 

s 

/ = ^/^V _i , /«£K for i = l,...,s. (4.4) 

Then we obtain 

s s s 

o = {hj} = y^ Us - * w - E ifijy"* = E f^ s - uyy 3 -'. 

i=0 i=0 i=0 

Hence, fi = for 2i 7^ s, and if for even s = 2r, f r 7^ 0, then / = f r (xy) r . This implies 
that every homogeneous, and thus arbitrary polynomial first integral / G K[x, y] of H , is 
an element of K[h], where h = xy. 

Now, assume that / is a rational first integral of H. Then we can write f = P/Q 
where P and Q are relatively prime polynomials in K[ar, y]. Hence we have 

= {H, P/Q} = (Q{H, P} - P{H, Q}) , 

so Q{H, P} = P{H, Q}. As P and Q are relatively prime, this implies that 

{H,P} = 1 P and {H,Q} = 7 Q. (4.5) 

for a certain 7 G K[x, y\. Comparing degrees of both sides in the above equalities, we 
deduce that 7 G K. If 7 = 0, then P and Q are polynomial first integrals of H, so in this 
case we have that / G M,(h). 

We show that the case 7 / is impossible. Let us assume that 7 / 0. It is easy to see 
that if P G K[x, y] satisfies equation 

{H,P}= 7 P, (4.6) 

then its every homogeneous component also satisfies this equation. Thus let us assume 
that P is homogeneous of degree s. If we write 

s 

P = Y^ PiX l y s '\ Pi G X for i = 1, . . . , s, (4.7) 

i=0 
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then equation (|4.6j> leads to the following equality 

s 

^P i (s-2i- 7 )xV~ i = 0. (4.8) 

i=0 

This implies that if coefficient Pi 7^ 0, then 7 = s — 2i and P = PiX l y s ~ l . Thus, every 
homogeneous solution of (|4.6|) is a monomial of the form PiX l y' l+1 , where 7 is a non-zero 
integer and i is a non- negative integer such that i + 7 > 0. Thus a non- homogeneous 
solution of (|4.6p is a finite sum 

i+7>0 

But Q satisfies the same equation (I4.6D . so we also have 

i+7>0 

If 7 > 0, then P and Q are not relatively prime, because they have a common factor y 7 . 
On the other hand, if 7 < 0, then they are not relatively prime, either, because they have 
a common factor x. We have a contradiction and this finishes the proof. □ 

Let us denote by 5 the Lie algebra which is the direct sum of n copies of sp(2, C) 

5 :=sp(2,C)e---©sp(2,C), (4.9) 

* v ' 

n— times 

and by 7Tj : s — > sp(2, C), the projection onto the i-th component of 5, for i = 1, . . . , n. 

Lemma 4.2. Let q be a Lie subalgebra of 5, and Qi = vrj(g) for i = l,...,n. Assume 
that q has 2n — 2 rational algebraically independent first integrals. If for some 1 < j < n, 
dim c Qj > 1, then Qi = {0}, for i / j. 

Proof. We consider elements of Lie algebra s as degree two homogeneous polynomials of 
2n variables (x,y). Because 5 is a direct sum of the form ()4.9p . we have 

rt 

s ~ 0C 2 [xi,2/i], 
i=i 

and an element H £ 5 admits a unique representation 

n 

H = '^2Hi, Hi G C 2 [xi,yi], for i = l,...,n. 

i=l 

Without loss of generality, we prove our lemma for j = n. Let us assume that dime n > 1, 
and let / G C(x,y) \ C be a first integral of q. Then from the proofs of Propositions 14.21 
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and 14.31 it follows that / does not depend on (x n ,y n ). Thus, if g has 2n — 2 algebraically 
independent rational first integrals fi, then 

fi eC(xi,...,x n _i,yi,...,y„_i)\C, i = 1, . . . , 2n - 2. (4.10) 

Assume that there exists < i < n such that Qi ^ 0. Then there exists a non-zero 
H G C2 [iCi,2/i] such that 

{/jiH}= !^_^ = „, for i= ,... 2 „_ 2 . 

<9x; dyi dyi dxi 

Hence, the linear Hamiltonian vector field Xh considered as a vector field in C 2n_2 with 
coordinates (x\, . . . ,x n -\,y\, . . . ,y n -i) has 2n — 2 first integrals (|4.1U|) and so Xh = 0. 
Thus we have a contradiction because for a non-zero H G C.2[xi,yi], Xh 7^ 0. □ 

Lemma 4.3. Let q be an Abelian Lie subalgebra of s, and Qi = vr.j(g) for i = 1, . . . ,n. 
Assume that q has 2n — 2 rational algebraically independent first integrals. Then dime < 
2. 

Proof. As g is Abelian, is Abelian, and dime Qi < 1 for i = 1, . . . , n. Let dime g = s. 
Then among Qi, exactly s have dimension 1, and the remaining are zero-dimensional. We 
can assume that dime Qi = 1 for i = 1, . . . , s, and dime Qi = for i = s + 1, . . . , n. 

Let / £ C(x) 3 , then, in particular / is a first integral of gi. We can consider / 
as an element of K(xi, y\), where K = C(x2, . . . , x n , yi, ■ ■ ■ > Vn)- Then, by Lemma 14.11 
there exists hi £ C(xi,yx) such that / G K(/ii) = C(hi, X2, • • • , x n , yi, ■ ■ ■ , y n )- Repeating 
successively the above reasoning for g^, for i = 2, . . . , s, we obtain that / is an element 
of R := C(hi, . . . ,h s ,x s+ i, . . . ,x n ,y s+ i, . . . ,y n ), where hi G C(xi,yi) \ C, for i = 1, . . . , s. 
Thus a first integral of g depends on r = s + 2{n — s) = 2n — s variables. This implies 
that every set of p > r elements from R is algebraically dependent. Hence, if /1, . . . , f p 
are algebraically independent elements of R, (i.e. algebraically independent first integrals 
of g), then p < r. So, from assumption of our lemma p = 2n — 2 < 2n — s, and thus 
s < 2. □ 



4.2 Proofs of Theorems D and [QI 

Under assumption of both theorems, the differential Galois group of the variational equa- 
tions is a direct product (|3.2I) . Thus we also have 

9° = 9{k, Ai)° x • • • x 9(fc, A n )° C Sp(2n, C). (4.11) 

Moreover, the Lie algebra g is a direct sum 

n 

= 00i, giGsp(2,C). (4.12) 
i=i 

We have already remarked that the identity components of the differential Galois groups 
of equations (|3.ip . (|3.3p and ()3.8[) are the same. Hence, conclusions of Lemma 13.11 and 
Lemma 13.21 are valid for S(&, A)°. 
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We prove Theorem 11.21 by contradiction. Let us assume that S° is not Abelian. Then 
there exists 1 < i < n such that 9(k, Aj)° is not Abelian. We have two possibilities. Either 
S(Aj)° is solvable and in this case it conjugates to the whole triangular group of Sp(2, C), 
or it is not solvable and then it is Sp(2,C), see Proposition 2.2 in [24j |. In both cases 
dimension of Lie algebra Qi of 9{k,Xi)° is at least 2. Now, by Lemma l4~2j for r / i, Lie 
algebras g r of groups S(k,X r ) are zero dimensional. 

We know that among Xj one, let us say A s , s € {1, . . . , n}, equals k — 1, thus it belongs 
to the item 1 of the table (j!.13f) . For A = A s the identity component G(k, X)° of differential 
Galois group of hyper geometric equation (|3.3p is Abelian but its Lie algebra has dimension 
one, see Lemma [3.21 But G(k,X)° = S(fe, A)°, so Lie algebra q s has dimension one. It 
follows that s 7^ %. But it is impossible because we have shown that for each k ^ i, 
k E {1, . . . , n}, Lie algebras Qk are zero dimensional. A contradiction finishes the proof of 
Theorem 11.21 

Now, we prove Theorem ll.3l From Theorem (II. 2D it follows that S° is Abelian and thus 
S(&, Aj)° are Abelian for i = 1, . . . , n. As 9(&, Aj)° = G(k, X{)°, we can apply Lemma [37T1 
and we obtain immediately the first statement of Theorem [L3j The last statement of this 
theorem follows from Lemma 14.31 and Corollary 13.11 



5 Global analysis 

The aim of this section is to present a pro of of Theorem 11.41 The central role in our 
considerations plays a theorem proved in ll| , see also 10(] . This theorem is related to the 
Kovalevskaya and Painleve analysis. To formulate it we recall basic facts. 
Let 

d_ 

dxi 



e(x) :=^Xi — , x = (xi, . . . ,x n ), (5.1) 



be the Euler vector field. A polynomial vector field v given by 

n d 

v(x) := } Vi(x) — , Vi£C[x], i = l,...,n, (5.2) 

r-f oxi 

i=i 

is called homogeneous of degree p, iff L e (v) = pv, where L e denotes the Lie derivative 
along the Euler field. Notice that a polynomial vector field v is homogeneous of degree p 
iff Vi are homogeneous and 

degVi=p+l for i = l,...,n. 
With vector field (|5.2p we associate the following polynomial differential equation 

^x = v(x), x£C n . (5.3) 

In this context, abusing notation, we consider v as an element of C[a;] n . We say that 
system (j5.3j) is homogeneous of degree p iff v as a vector field is homogeneous of degree p. 
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Definition 4. We say that a direction d £ C n is a Darboux point of a homogeneous 
vector field v, or a Darboux point of system (|5.3p . iff e(d) A u(d) = and v(d) 7^ 0. 

Thus vector field v and the Euler field are parallel at a Darboux point. This implies 
that there exists a unique p £ C* such that vector field pv — e has an equilibrium point 
at d. To find all Darboux points of a given vector field v, it is enough to find all non-zero 
solutions of non-linear equations 

v(d) = d, (5.4) 

and distinguish among them those giving different directions. We always assume that a 
Darboux point d is normalized in such a way that it satisfies (15. 4p . Notice also that a 
Darboux point d = (di, . . . , <i n ) can be considered as a point of the projective space CP n ~ 1 
with homogeneous coordinates [d\ : ■ ■ ■ : d n ] . 

The set of all Darboux points of a vector field v we denote by D v . It can be empty, 



finite or infinite. One can show, see Proposition 4 in [111 ], that if all Darboux points of a 
homogeneous vector field v of degree p are isolated, then their number is not greater than 
d(n,p) := ({p + l) ra — l)/p, and a generic v has exactly this number of Darboux points. 

Let d be a Darboux point of a homogeneous vector field (15. 2D . The Kovalevskaya 
matrix K(d) at point d is defined as 

#(d) := «'(d) - £, (5.5) 

and its eigenvalues are called the Kovalevskaya exponents. 



Remark 5.1 In the above definition we assumed that Darboux point d satishes (jBT 
Without this assumption the Kovalevskaya matrix at Darboux point d is dehned as the 
matrix of linearization of vector field pv — e at equilibrium point d, i.e., 

K(d) := pv'{d) - E. 

We have the following fact. 

Proposition 5.1. Let d £ C n be a Darboux point of homogeneous vector field v of degree 
p. Then d is an eigenvector of the Kovalevskaya matrix K(d), and the corresponding 
eigenvalue is p. 

An easy proof of the above proposition we leave to the reader. 

Let Aj = Aj(d) for i = 1, . . . ,n denote the Kovalevskaya exponents at a Darboux point 
d £ T> v . By Proposition [5J3 one of them is p. We always assume that A n = p, and denote 
by A = A(d) = (Ai, . . . , A n _i) £ C n_1 . Let r, for < i < n — 1, denote the elementary 
symmetric polynomial in (n — 1) variables of degree i. In paper [ll| it was shown that 
the values of A(d) taken at different Darboux points are not arbitrary, i.e., there exist 
universal relations among all A(d), d £ T> v . 

Theorem 5.1. Assume that a homogeneous polynomial vector field v of degree p > has 
d(n,p) isolated Darboux points, and let S be a symmetric homogeneous polynomial inn—1 
variables of degree lower than n. Then, the number 

^-i(A(d)) 
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depends only on a choice of S, the dimension n and the degree p. 



In other words, R = R(S,n,p) does not depend on a specific choice of v, provided 
that v has the maximal number of Darboux points d(n,p). For applications of the above 
theorem we have to know values R(S,n,p) for an arbitrary n, p, and for a chosen set of 
independent symmetric polynomials S of degree not greater than n — 1. The standard 
choice for S are elementary symmetric polynomials Tj, or alternatively t[, for < i < n— 1. 
To calculate R(S, n,p) it is enough to find a homogeneous system defined for an arbitrary 
n > 2 and p > possessing the maximal number of Darboux points, and such that 
one can determine the Kovalevskaya exponents easily for all its Darboux points. These 



requirements are satisfied be a ra-dimensional generalization of the Jouanolou system 18|] 



ii = xf+i, 1 < i < n, x n+1 = Xi. (5.7) 

For this system, the Kovalevskaya exponents at all Darboux points are the same and are 
given by 

Aj = (p + l)e n_i - 1, 1 < i < n, (5.8) 



where e is a primitive n-th root of the unity. In 31] we have shown that 



r r (A) = (-ir^( n r ^ \p + l)\ (5.9) 

i=Q ^ ' 

for < r < n — 1. Using this fact one can show that 

R(^,n,p) = (-ir +r - l (n + py, (5.10) 

and 

R(r r , n,p) = (-1)™— 1 £ ( U l ~) (p + 1)* , (5.11) 

for < r < n — 1. 

Now we return to the Newton equations (|1.5|) . We associate with them the following 
auxiliary system 

— x = F(x), x£C n . (5.12) 

As degFj = k — 1 for i = 1, . . . n, F considered as a vector field has degree p = k — 2. 
A Darboux point of force F is a Darboux point of the auxiliary system. Thus, if d is a 
Darboux point of (I5.12p . then among the Kovalevskaya exponents calculated at d, one is 
k — 2. We denote the remaining ones by A(d) = (Ai(d), . . . , A n _i(d)). Assuming that 
the auxiliary system has d(n,k — 2) = D(n,k) isolated Darboux points we can apply 
Theorem 15. 11 Then, taking into account relations (|5.10p and (|5.1ip . we obtain the thesis 
of Theorem 11.41 and this finishes its proof. 

One of the most important consequences of Theorem 11.41 is the following. Let J-(n, k) 
denote a set of forces with exactly L = D(n, k) Darboux points dj, j = 1, . . . , L. For each 
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F G J~(n, k) we have k = (n — 1)L nontrivial Kovalevskaya exponents. We denote their 
collection by K(F) G C K , i.e., 



K(F) := (A(di), . . . ,A(d L )). 



If for fC{F) given above, all (k, Ajj), Xij = Aj(dj) + 1 for i = 1, . . . , n— 1 and j = 1, . . . , L 
belong to table (j!.13[) in Theorem 11.31 then we write K,{F) G J(n,k). In other words, if 
for F G J-(n, k) the corresponding Newton's equations are integrable in the Jacobi sense, 
then JC(F) G J(n,k), and if JC(F) G J{n,k), then the Newton equations satisfy the 
necessary conditions for the integrability in the Jacobi sense of Theorem 11.31 We show 
that sets J(n, k) are finite. 

Theorem 5.2. For n, k > 2 set J{n, k) is finite. 

Proof. Relation (|1.15|) in Theorem 11.41 for r = n — 2 and p = k — 2 reads 



By definition of J(n, k), quantities Xij = Aj(dj) + 1 for i = 1, . . . , n — 1 and j = 1, . . . , L 
belong to sets listed in table (|1.13p in Theorem 1 1.31 Because of this, we can apply directly 
Lemma B.l from [22| and this finishes the proof. See also the proof of Theorem 4 in 



6 Two degrees of freedom 

In this section we investigate in detail the case of the Newton equations with two degrees 
of freedom. The main aim is to show that in this case we can obtain more general results 
than those given by Theorem 11.41 Moreover, we investigate also the case of forces which 
admit infinitely many Darboux points, as well as cases when forces do not have Darboux 
points. 

Before we start our analysis of the system with two degrees of freedom, we make two 
general remark. If a force admits Darboux points d%, ■ ■ ■ , di, 1 < I < n, then we can choose 
I coordinate axes along them. In a such choosen frame in such all but j-th coordinates of 
dj are zero. 

Newton's system of the form (jl.5p contains as a subclass natural Hamiltonian systems. 
Thus it is convenient to have a simple criterion to distinguish such systems. In Appendix 
A we present two methods that enable to check the equivalence of the analyzed Newton 
system to a Hamiltonian one. 

In the rest of this section we assume that F = [F\ , F<i ) , and we put 




(5.13) 



□ 



fc-i 




(6-1) 
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6.1 Darboux Points 



Darboux points of a force F are the equilibria of the vector field G(q) = F(q) — q. Thus, 
it is convenient to introduce an auxiliary system 

qi = -qi + F 1 , q2 = -q2 + F 2 . (6.2) 

However, we must remember that different equilibria of the above system can correspond 
to the same Darboux point. To cope with this problem we introduce coordinates x = q% 
and z = q 2 /qi and define 

fi(z):=Fi(l,z) for i = l,2. (6.3) 
Then the auxiliary system (|6.2|) in these coordinates takes the form 

x = -x + x k ~ l h(z), z = x k ~ 2 g{z), (6.4) 

where 

h(z)=f 1 (z), g(z) = f 2 {z)-zh{z). (6.5) 

Notice that degg < k. As it has been mentioned, a Darboux point is a point in a 
projective space, and here it is CP 1 with \q\ : q 2 ] as homogeneous coordinates. Hence, z 
is a coordinate on the affine part of CP 1 . An equilibrium of (|6.4|) is a solution of 

x k - 2 h(z) = 1, x k - 2 g(z) = 0, x^O, (6.6) 

and thus the Darboux points are roots of polynomial g(z) such that h(z) does not vanish 
on them. This immediately implies that either the number of Darboux points is infinite 
when g = 0, or it is finite and not greater than k. We say that a Darboux point is multiple 
iff its z coordinate is a multiple root of g(z). 

Remark 6.1 As it has already been mentioned, not all roots of g(z) are Darboux points. 
If g(z±) = h(z+) = 0, then z± is not a Darboux point. Equations g(z±) = h(z*) = imply 
that fi(z+) = f 2 {zi,) = 0, i.e., F\ and F 2 have a common factor P = q 2 — z+q\. 

If [0 : 1] € CP 1 is a Darboux point, then we say that it is located at the infinity. To 
investigate its neighborhood we use coordinates y = q 2 and C = Qi/q2- We define 

fi(Q :=F 4 (C,1), for z = 1,2, (6.7) 

and 

g(0 ■= /i(C) - CMC), HO ■= MO- (6.8) 

Then a Darboux point is a root of polynomial g, such that it is not a root of h. If 
5(0) = /l(0) = and simultaneously h(0) = f 2 (0) / 0, then the infinity is a Darboux 
point. Notice that 

/i(C) = C*~Vi Q) for i = l,2. (6.9) 
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In many cases, it is convenient to assume that all Darboux points are located in the 
affine part of CP 1 , i.e., that there is no Darboux point at the infinity. If degg = k, then 
surely the infinity is not a Darboux point. Hence, let us analyze whether we can assume 
that for a given force the corresponding polynomial g is of degree k. Prom the definition 
of polynomial g, see (I6.5|) . it follows that degg = k if and only if deg/i = k — 1. We show 
the following. 

Proposition 6.1. Let F = (iq,^), F\,F 2 G C[(/i,Q2] be a homogeneous force and 
degFi = k — 1 for i = 1,2. Assume that for each A G GL(2,C) for the equivalent 
force Fa = {Ff, F^) we have deg ?2 F-f 4 < k - 1. Then 

G := qi F 2 - q 2 F 1 = 0. (6.10) 

Proof. One can directly check that if F satisfies condition (|6.10j) . then every force equiv- 
alent to F also satisfies this condition. 

First let us notice that condition deg 92 iq < k — 1 is equivalent to deg fx < k — 1, and 
this is equivalent to /i(0) = 0. In fact, if 



fi(z) = Y,fh\ fi^O and Kk-1, 



then 



and hence /i(0) = 0. 
Let 



Then 



i=0 



/i(0 = c fc -7i 



A 



a (3 
7 5 



(7) 

V s / i=0 



A = aS + /fy / 0. 



Fi (Qi,Q2) = 2£ [8Fi(aqi + /3q 2 ,jqi + Sq 2 ) 



and hence 



where 



(3F 2 (aqi + /?<?2,7<?i + 5q 2 )\ 



(6.11) 



/i 4 (C) = i(7C + ^) 



k-l 



r(0 



Sfi(T(0)-Ph(r(C)) 
aQ + (3 



7 C + <5 

Assume that 5^0 and denote x = (3/5. Then we obtain 



ff(0) = -[f 1 (x)-xf 2 (x) 

By assumption /^(O) = for an arbitrary A, so 

fi(x) - xf 2 (x) = 0, 



(6.12) 



for an arbitrary x G C. Putting x = qi/q 2 in the above equality an using the definition of 
fi and f 2 , see (|6.7p . we obtain condition (|6. 10j) . □ 
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Remark 6.2 Condition (16. 101) means that force F = (Fi, F2) admits infinitely many Dar- 
boux points. 



6.2 Generic case 



Now we determine the Kovalevskaya exponents at a given Darboux point. The Jacobian 
matrix of the auxiliary system (|6.4p has the form 



J 



-1 + (k - l)x k ~ 2 h(z 



„k-lu' 



h'{z) 



(k — 2)x k 3 g(z) x k 2 g'(z) 
For a Darboux point (x, z) satisfying (|6.6[) this matrix becomes 



(6.13) 



J 



k — 2 xh'(z)/h(z) 
g'(z)/h(z) 



Hence the non-trivial Kovalevskaya exponent at a Darboux point with coordinate z is 
given by 

g'(z) m - h(z) - zf[(z) 



A(z) 



(6.14) 



Kz) h(z) 

On the projective line CP 1 , where Darboux points are located, we define a meromorphic 
differential form uj which in coordinate z is given by 



h(z) 
9(z) 



dz, 



and in coordinate <^ = l/z in the neighborhood of infinity by 

h(i/c) dc 



9(1/0 C 



2 ' 



(6.15) 



(6.16) 



Notice that if z+ is a Darboux point, then it is a pole of uj, and the residue of uj at this 
point is 



res (a;, z+ 



1 



(6.17) 



9'{z*) k{z*y 

provided that g and h are relatively prime and z* is a simple root of g. 

For a given force F = (Fi,^) with finitely many Darboux points, we can always 
assume that all of them are located in the affine part of CP 1 i.e. for all Darboux points 
q± ^ 0. This can be achieved by a proper choice of coordinates. Under this assumption 
we have Fi(0,q2) 7^ 0, and hence 7^ 0, see (|6.1[) . We keep this assumption in our 

further considerations. 

To formulate a generalization of Theorem ll.4l we have to distinguish a class of forces for 
which this is possible. Let Ai denote a set of forces F = (F\,F2) satisfying the following 
conditions 



A.l F admits 1 < I < k simple Darboux points; 



23 



A. 2 if F\{q) and G{q) := 92^1 (q) — Qi^iq) have a common linear factor P = aq\ + /3q2, 
M + \0 I 7^ 0) then the multiplicity of P in G is not greater than the multiplicity of 
P in Fx. 

Theorem 6.1. Assume that F e *4; ; and let A, is i/ie non-trivial Kovalevskaya exponent 
of F at i-th Darboux point for i = 1, . . . , I. Then 

Ex: = - L (fU8) 

1=1 

Proof. We can assume that deg<?(z) = k and deg/t(^) = k — 1. Moreover, from the 
definitions of h and g, see (|6.5p . it follows that if az k ~ l is the highest order term of h, 
then — az k is the highest order term of g. 

Although for I = k the statement of theorem coincides with Theorem 11.41 it is instruc- 
tive to give first an independent proof for this case. Thus, let k = I. Then g and h are 
relatively prime and all k roots Zi of g are simple. We have 

To calculate the residue of a; at infinity, i.e., res(a> oo ,0), we notice that for polynomials 
h(C) = C, k - l h{l/C) and g(Q = C, k g{l/0 we have h(0) = a and g(0) = -a. Thus 

h(i/o _ Mo _ 1 , 



C 2 5(i/C) CKC) C 

where dots denote higher order terms. This shows that res(a; 00 , 0) = 1 and the application 
of the global residue theorem, see e.g. [jj, finishes the proof for I = k. 

For 1 < I < k the assumption that F £ Ai implies that g(z) has I simple roots 
zi, . . . , Zi, such that h{z{) 7^ for i = 1, . . . , I. Moreover, if s {z\, . . . ,Zi} is a root of 
g(z) with the multiplicity j, then s is also a root of h{z) with the multiplicity greater or 
equal to j. All these conditions imply that we can write 

g (z) = -ap(z)g(z), h(z) = ap(z)h(z), 

where p{z), g(z) and h(z) are monic polynomials; p is of degree k — I, a G C* and 

I 

K z ) = J[( z - z i)- 



i=l 



Hence 



h(z) 

g( z ) 



and then 

1 h(zi) p(zi)h(zi) h(zi) 



res((j,Zj), (6.19) 



K g'( z i) P'(zi)g(zi) +p(zi)g'(zi) g'(zi 
for i = 1, . . . , I. Moreover, it is easy to see that res(u;, 00) = 1. This finishes the proof. □ 
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Remark 6.3 If F satisfies assumption A.l and does not satisfy A. 2 we still can do some- 
thing. In such situation we have 

s I 

g(z) = -ap(z)g(z), p{z) = \[{z - Si ) n \ ~g{z) = \[(z - z t ), (6.20) 

and 

s 

h{z) = aq(z)h(z), q(z) = [](z - Sl ) m > , (6.21) 

i=i 

where 

s s 

I + rii = k, deg h + mi = k — 1, 

i=l i=l 

Si {zi, . . . , z{\ and h(si) / for i = 1, . . . , s. Now 

oj = - .. , n . r dz, (6.22) 

and the global residue theorem gives 

l ^ s 

EAT + ^C res(w ' Sj) = " L (6 - 23) 

i=l 1 j=\ 

We notice that if rii > rrii, then, generally, residue res(w,Sj) / and it depends on 
coefficients of polynomials F\ and F2. 



6.3 Non-maximal number of Darboux points and multiple Darboux 
points 

If a force F = {F\,F2) does not have the maximal number of Darboux points, then 
polynomials F\ and F2 have a common factor P £ C[q] \C, i.e., Fi = PFi, where Fi £ C[q] 
for i = 1,2. In such cases it is possible to find a particular solution which does not 
correspond to a Darboux point. For example, if P = q2, then qi(t) = jt, q2(t) = 
is a particular solution. Amazingly, we can use such simple solutions and methods of 
differential Galois theory for proving the non-integrability. First we show a very simple 
fact. 

Proposition 6.2. Let F(q) = (a\qi + ot2q2){Fi{q),F2{q)) where \ot\\ + |ck2 | / 0, and 
Fx, F 2 e C[q]. Then there exists A £ GL(2,C), such that AF(A~ l q) = q 2 (F 1 (q), F 2 {q)), 
where Fi £ C[q] for i = 1,2. 

Proof. Matrix A given by 

det A = I ct2 1 2 + 1 0:2 1 2 , 

is non-singular and satisfies the desired condition. Here asterisk denotes the complex 
conjugation. □ 



A 



a 2 
ai 



OL2 
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Hence, if a force F = {F\,F%) where deg-Fj = k — 1, has less than k Darboux points, 
then we can always assume that F\ and F2 have a common factor qi- We show that under 
a certain genericity condition the Newton equations with such force are not integrable in 
the Jacobi sense. 

Theorem 6.2. Assume that F±, F2 € C[qr] are homogeneous polynomials of degree k — 2. 
If deg gi F2 = k — 2, then system 

9i= Pi, Pi = (91,92), 92=P2, P2 = -92-^2(91,92), (6.24) 

does noi admit two functionally independent polynomial first integrals. 
Proof. System (16.24)) possesses the following particular solution 

q 1 (t)= 1 t + (3, pi(t)=7, 92(0 =P2(t) =0. 
Variational equations along this solution take the form 



We have 



(Qi\ 




f 





1 









(Qi\ 


Pi 






-Fi(qi(t) 













Pi 


Q2 















1 




Q2 


w 












-F 2 (qi(t),0) 


0) 






^1(91,0) 






and 


F 2 (qi(t),0) 


= a29i 2 , 



(6.25) 



where «2 / 0, as, by assumption deg gi ^2(91, 0) = k—2. We choose constants of integration 
(3 = and 7 in such a way that a27 fc_2 = — 1- Then the variational equations take the 
form 

Qi=Pi, P 1 =at m Q 2 , Q2 = P2, P2=t m Q 2 , «eC, m = k-2 > 0. (6.26) 

First we show that the differential Galois group G of the above variational equations is 
a direct product {E} x SL(2,C), where E is the 2x2 identity matrix. To this end, we 
notice that Q2 satisfies the following equation 

Q2 = t m Q 2 . (6.27) 

Let K D C(t) be the Picard-Vessiot extention for this equation, i.e., K is the smallest 
differential field containing two linearly independent solutions of (|6.27p . Then all solutions 
of (|6.26p belong to K 4 . It is clear for a = 0. If a 7^ 0, then Qi — aQ2 = 0, and so Qi(t) = 
&Q2(t) + at + b, for some a, b G C. So, as we claimed (Qi(t), Qi(t), Q2{t), Q2(t)) S K 4 . 

Now we can determine the differential Galois group of variational equations (|6.26[) . For 
this purpose we have to know the differential Galois G group of equation (|6.27j) . In [22[] it 
was shown that G = SL(2, C) for an arbitrary m > 0. Let w\ and u>2 denote two linearly 
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independent solutions of (|6.27p . Then the fundamental matrix of (16.260 has the following 
form 

1 t a.w\ aw<i 
1 aw\ aw2 

W\ W2 
W\ W2 



X 



(6.28) 



Now, if for g G G we have 



g(wi) = cji^wi + a 2 ,iW 2 for i = 1, 2, 



(6.29) 



then £ = [crjj] G SL(2,C), and 



g{X) = XS, S 



E 0' 

s 



(6.30) 



where is the 2x2 identity matrix. Hence the differential Galois group G of system (I6.26P 
is a direct product {E} x SL(2,C). 

Let us notice that we cannot finish our proof at this point invoking Lemma [2.1l In fact, 
the differential Galois group G has two independent invariants Q\ and P\. Thus we proceed 
in the following way. From the proof of Lemma l4,2l it follows that if R £ C(Q\,Pi, Q2, P2) 
is an invariant of G, then R G C(Q\, P±). This implies that if / G C(t)(Qi, P±, Q2, P2) is 
a first integral of variational equations (|6,26[) . then / G C(t)(Qx, -Pi). 

Let us assume that P G C[q\,px, (72, P2] \ C is a first integral of (|6.24p . We can write it 
in the following form 

P= Yl P i,o(^Pl)<AA- (6-31) 

0<i+j<m 

The leading term p G C[t](Qi,Pi, Q2, P2) of P is the lowest degree term of 

p(rrt + Q 1 , 1 + p 1 ,Q 2 ,p 2 )= Yl PiA^ + Qi^ + p i)Qi p i, ( 6 - 32 ) 

0<i+j<m 

and it is a first integral of variational equations (I6.26p . We have already shown that the 
integral p must belong to C[t](Qi, Pi). Hence p is the lowest order term of Pq o(7*+Qi ; 7+ 
Pi). 

We show that Po t o G C[p{\. Inserting expansion (|6.3ip into equation 



dP dP dP dP 

Pit. ^P2j. q2Fi(qi,q 2 )j, ?2-^2 (91,92)^— = 0, 

Oq\ oq 2 dpi dp 2 



(6.33) 



obtained from the condition P = 0, and taking terms of degree zero with respect to (92,^2) 
we obtain 

dP ,o n 
Pi-z — = 0, 



dqi 



and hence Pq,o ^ as we claimed. 
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Now, assume that system (I6.24D admits two independent polynomial first integrals P 
and Q. We can assume that their leading terms p,q £ C[t](Qi, Pi, Q2, P2) are independent. 
But we know that p,q S C[t](Qi, P\), and thus, p and q are the lowest degree terms of 
Pofiilt + Qii 7 + A) and Qo,o(7^ + Qi-, 7 + Pi)- However, we have already shown that Po,o 
and Qo,o depend only on pi, so p and q depend only on Pi, and thus they are functionally 
dependent. A contradiction finishes the proof. □ 

Remark 6.4 If in the above lemma we assume that deg qi F\ = deg qi F% = k — 2, then a in 
variational equations (|6.26p is not zero, and we can modify the final step of the proof. Let 
P £ C(t)(Qi, Pi) be a first integral of variational equations (16.260 . Then the condition 
P = yields 

dP n dP m ^ dP n 

dP 

and thus — — = 0. This implies that the variational equations admit at most one rational 
0P1 

first integral, and thus the considered system (|6.24p is not integrable in the Jacobi sense 
with two rational first integrals. 

Remark 6.5 Let us assume that F{q) = (aiqi + 02^2) (Fi(q), F 2 {q)) ■ Without loss of 
generality we can assume that 02 7^ 0. If 02 = 0, we just renumber the coordinates. By 
Proposition \6.2\ there exists a non-singular matrix A such that 

F A (q) := AF{A- l q) = q 2 (Fi(q), F 2 (q)). 

Assumption deg qi F 2 = k — 2 in Theorem \6.2\ means that ^2(51, 0) 7^ 0. We have 

■^2(91,92) = aiFi(a 2 qi + 0^92, -«i9i + "292) + a2-^2(a29i + «i92, -aiqi + a%q 2 ) , (6.34) 
and thus 

F 2 (qi,0) =af- x q\~ 2 aiFi(a 2 , -a x ) + a 2 F 2 (a 2 , -«i) . (6.35) 
The above shows that if the polynomial 

G(q) := qiF 2 (q) - q 2 F~i (q), 

is divisible by (ai9i + 0:2(72), then F 2 (qi,Q) = 0. On the other hand, if F 2 (qi,Q) = 0, then 
from (|6.35p it follows that 

^ ^ Oti 
g(a) := i^fl, a) — aFi(l,a) = /2(a) — afi(a) =0, a = . (6.36) 

OL2 

Thus a is a root of polynomial g G C[z], and this is equivalent to the fact that (aiQ , i + a292) 
is a factor of G(q). 

It can happen that the considered force does not have Darboux points. Such forces 
are characterized in the following theorem. 
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Theorem 6.3. Assume that a force F = (Fi,^), degFj = k—1, does not have a Darboux 
point. Then 

k-l 



Fi = Y\(q2- mi), 



i=l 



F 2 = i 
Qi 



q S F 1 -l[(q 2 -z i qi) n ' 



i=l 



(6.37) 



where 1 < p < k—1, Zi 6 C, and z\, . . . , z p are pairwise different; n« G N and n\ + - ■ -+n p = 
k. Moreover, if there exist 1 < j < p such that rij = 1, then the Newton equations with 
force F do not admit two functionally independent polynomial first integrals. 

Let us note that in spite of the fraction in the expression on F2 the force F is polyno- 
mial, see the proof below. 

Proof. Without loss of generality we can assume that deg^ F\ = k — 1. Then for polyno- 
mials h, g £ C[z] defined by (|6.5f) we have deg g = k and deg h = k — 1. Assumption that 
F have no Darboux point implies that all roots of g are also roots of h. So g has at least 
one multiple root and we can write it in the following form 

p p 
g(z) = a Y[(z - Zi) n \ rtj G N, = k, a £ C*, 



i=i 



i=i 



where 1 < p < k — 1 and z±, . . . , z p are pairwise different roots of h. Let z p+ i, . . . , z^-i be 
the remaining (not necessarily different) roots of h. Then we can write 



k-l 



i=l 



h(z)=pn 

Thus f\(z) = h(z), and 

f 2 (z) = g(z) + zh(z) = a ]J(z - Zi ) n * + (3z\{{ 



k-l 



Z - Zi 



i=l 



i=l 



But deg/2(z) < k, thus a = —f3 and we can put f3 = 1. Now, as Fi(qi, 52) = q\~ l fi (92/91) 
for % = 1, 2, we obtain (|6.37|) and this finishes the first statement of our theorem. 
To proof the second statement let us assume that n\ = 1. Then we can write 



Fi = - ziqi)Fi for i = 1,2, 



where 



fc-i 

Fi = JJfe - mi) 

i=2 



F 2 = - 
1i 



Q2F! - Yl(q 2 - mi) ni 



i=2 



Let us calculate polynomial g(z) as in Remark 16.51 We have 



tj(z) = G(\,z) = F 2 (l,z)-zF 1 (\,z) 



i=2 



Z - Zi 



(6.38) 
(6.39) 

(6.40) 



so g{z\) 7^ 0. Thus, as we explained in Remark 16.51 one can apply Theorem 16.21 to the 
considered case and this finishes the proof. □ 
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Let us notice that if the assumptions of the above theorem are satisfied, then the 
considered Newton's system still can possess a single first integral. Let 

fc-i 

■Pl(?l>?2) = IJ(?2 - mi) and F 2 (qi,q 2 ) = z j F 1 (q 1 ,q 2 ), 

i=l 

where 1 < j < k — 1 and z\, . . . , z^-i are pairwise different. This force has the form (|6.37p 
with p = k — 1, Hi = 1 for i ^ j and rtj = 2. However, the respective Newton system 
possesses a first integral 

I = ZjPi -p 2 . 

If in (|6.37p all are greater than one, then the system can be integrable. We give such 
examples in the next section. 

6.4 Infinitely many Darboux points 

If 

91-^2(91,92) = 92^1(91,92), (6.41) 

then force F = (F%,F 2 ) has infinitely many Darboux points. At each of these points the 
Kovalevskaya exponents are A\ = and A2 = k— 1.. Thus the Morales-Ramis theory does 
not give any obstructions for the integrability. 

It is easy to show that condition (|6.41|) is satisfied if and only if 

fe-2 fc-2 

i^E/W?- 1 -^ and F 2 = J2fi<lt 2 ~ i <fr 1 - (6-42) 

i=0 i=0 

The Newton system with a force satisfying (|6.41|) possesses one first integral 

I = HP2 ~ 92P1- 

For k = 3 and A; = 4 the Newton equations with the above form (|6.42p of forces are 
integrable with an additional first integral quadratic in the momenta 

h = KhPi + hV2? + 2(/o9i + /192) 3 , 

h = 2(/ p? + /1P1P2 + hp\) + {fml + /19192 + /29D 2 , 

respectively. But for k > 4, it seems that the Newton equations with forces given by (|6.42|) 
are not integrable in a generic case, nevertheless, they admit an additional first integral 
for some special values of /j. For example, if 

Fx = qiifoqi + fm) 1 " 2 and F 2 = q 2 (f oqi + f iq2 ) k ~ 2 , (6.43) 

then the additional first integral has the following form 

I = Hfopi + hV2? + 2(/ 9i + fm) k . (6.44) 

We do not know how to investigate systematically the integrability of the Newton equations 
in this case. It seems that we can apply only a direct method. 
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Lemma 6.1. Assume that force F = (F±, F2) admits infinitely many Darboux points and 
the corresponding Newton equations have a polynomial first integral which is a quadratic 
form with constant coefficients with respect to the momenta. Then these equations are 
equivalent to an integrable Hamiltonian system with the radial potential 

H =\{p\+pl) + \{ql+ql) k ^ 
for k even, or to the non- Hamiltonian integrable system 

Fi = qiqt 2 , F 2 = qt\ h = \pI + \qI 

for any k. 

Proof. Without loss of the generality one can postulate the following form of the first 
integral quadratic in the momenta 

k 

I = ap?+p2+J>g*-V& (6.45) 

i=0 

The case when the quadratic form of momenta is nonsingular corresponds (after rescaling) 
to a = 1 and the case when this form is singular corresponds to a = 0. Condition 1 = 
on solutions of (j6.42j) yields 

(k - i)bi - 2afi = 0, i = 0, . . . k - 2, b k ^ = 0, 
6i = 0, (i + 2)b l+2 - 2f t = 0, i = 0,...,k-2, 

and the detailed analysis of solutions of this linear system yields the desired result. □ 



7 Two degrees of freedom. Applications 

In this section we consider also Newton equations with two degrees of freedom but now we 
apply general theorems from the previous sections to general homogeneous forces of low 
degrees. The first aim is to show that our general results allow to perform the integrability 
analysis in a systematic and algorithmic way. Let us describe the main steps of this 
algorithm. 

We fix k > 2 and consider the Newton system with a homogeneous polynomial force 
F = (F\,F2) where degiq = degi^ = k — 1. The problem is to give a complete list of 
all F for which the corresponding Newton system is integrable in the Jacobi sense. As 
C-linear space Cfc_i[q] of homogeneous polynomials in two variables of degree k — 1 has 
the dimension k, the problem depends on 2k parameters. We can reduce the dimension 
of the problem taking into account the action of the four dimensional group GL(2,C). 

If for a given k a complex number A belongs to an item in table (|1.13|) of Theorem 1 1.'6\ 
then we write A S A4- 
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In the first step of the algorithm we have to find for 1 < I < k all solutions of equation 



1 1 



i 1 



22 



such that Aj = Aj + 1 G A4 for i= 1, . . . , I. The set of such solutions we denote as in 
by 1^,1- Let us remark here that for small k it can be determined easily with the help of 
a computer algebra program, however, for bigger A: it is a time demanding computational 
problem. 

In the next step for each solution (Ai, . .. ,Aj) obtained in the previous step we have 
to determine all non-equivalent forces F with exactly I simple Darboux points, and such 
that at i-th point the non-trivial Kovalevskaya exponent is Aj, for i = 1, . . . , I. For small 
k it can be done effectively, however, for bigger k the problem is hard, because we have to 
solve a system of nonlinear polynomial equations of high degrees. The forces selected in 
this way satisfy all necessary integrability conditions of Theorem 11.31 but of course not all 
of them are integrable. To decide which of them are integrable, we have to use a theorem 
giving the stronger necessary integrability conditions than those contained in Theorem ll.31 
However, as far as we know, there is no such a theorem. In Appendix B we present a 
possible approach to this problem. It is a dire ct g eneralization of the method based on 



higher order variational equations developed in [2J; |27|] for Hamilton equations. 

The above steps do not finish the procedure. We have to investigate separately all 
the cases when the force have lower than the maximal number of Darboux points with 
a multiple Darboux point or does not have such a point at all, because for those cases 
we cannot apply Theorem 11.31 but, nevertheless, we can use Theorem 16.21 Moreover, 
in all considered cases we have to distinguish all instances when the Newton system is 
Hamiltonian. 

Below, in a case when we cannot apply any theorem to check the integrability, we use 
several different algorithms for finding a polynomial first integrals of a given degree. We 
do not describe details of these computations and below we call such a search a direct 
method. 



7.1 Case k = 3 

Here we investigate the Newton equations with a homogeneous force F = (i 7 !,^) of 
degree two 

rp a ll 2 , , fl 13 2 jp °21 2 , , a 23 2 /~ -, \ 

F\ = ~Y<ll +0129192 + ~Y<l2i F 2 = -^-9l + ^229192 + ~Y<l2- I 7 - 1 ) 

Let us assume that F possesses at least one Darboux point. Then it can be placed at 
the infinity, so i*i(0, 92) = 0, and thus in (|7.ip we put 013 = 0. Additionally, by a linear 
change of variables 91 = aQ\ and q2 = bQi + CQ2 where ac 7^ one can obtain 022 = 0. 

As explained in Appendix A, for k = 3 the Newton equations are almost always 
Hamiltonian. After the above simplifications the Newton system is not Hamiltonian only 
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when 



a\2 = and 021023 / 0, or (7.2) 

a 2 \ = and 012(^23 - 012) 0. (7.3) 

Hamiltonian systems with homogeneous potentials of degree 3 were fully analyzed in 
[20I ]. so we consider here only non-Hamiltonian cases. At first, let us assume that condi- 
tions (|7.2p are satisfied. Then, as 021 7^ and 023 7^ 0, we can make a rescaling which 
gives 021 = — 1 and 023 = 1- Hence we have to consider the following one parameter family 
of Newton's equations 



q\ = -^aqj, a G C, 



1 2 1 2 



(7.4) 



This system has one first integral 



I = 3pi + aq(, (7.5) 



and for the complete integrability one additional first integral is necessary. 
For (|7.4p polynomials g{z) and h(z), see definitions (|6.5|) . have the forms 

SO) = \(z 2 + az ~ !)> h ( z ) = 2 a - 

Thus 



a ± V« 2 + 4 
*M = 2 . 

are two Darboux points provided a 7^ 0. If a 7^ ±2i, then z\ 7^ ^2- The non-trivial 
Kovalevskaya exponents at these points are 



A M = ±yi + ± (7.6) 

Now, if system (j7.4j) is integrable, then Ai and A2 belong to table (11.130 of Theorem 1 1.31 It 
is easy to notice that there are only two possibilities: either Ai = A2 = 0, or Ai = — A2 = 1. 
But for Aj given by (|7.6p the second case is impossible. The first case gives a = ±2i. 

Summarizing, we have to investigate equations (|7.4p only when a = ±2i and a = 0. 

Assume that a = 0. We introduce new coordinates 

Qi = qi + <?2, Q2 = qi - Q2- 
In these coordinates equations (|7.4j) read 

Qi = -^QiQ2, Q2 = —-^QiQ2- 



33 



Applying Theorem 16.21 to this system, we conclude that for a = equations (17. 4D are not 
integrable in the Jacobi sense with polynomial first integrals. 

If a = ±2i, then the considered force has a double Darboux point. At this Darboux 
point the non-trivial Kovalevskaya exponent is A = 0. Making the following change of 
variables 

^ 1 ^ a f f-7 i-r\ 

Qi = —rqi + 2<?2, Q2 = ^qi, (7.7) 

we transform system (I7.4D for the prescribed values of a into the following one 

Q 1 = -Q 1 (Q 1 + Q 2 ), Q 2 = -Ql (7.8) 

and now the double Darboux point is at infinity. Let us consider a phase curve r e , e ^ 0, 
corresponding to this Darboux point. It is given by 

e = \Q\ + \Ql Qi = Qi = 0. 

The particular solution corresponding to this phase curve is given by Q 2 (t) = —6p(t, g2,gs), 
where invariants of the Weierstrass function are g 2 = and g^ = — e/18. Variational equa- 
tions along this particular solution have the following form 

£i = 6p(t,92,93)£i, 6 = 12p(i,52,53)6, (7.9) 

i.e., they are a direct product of two Lame equation of the Lame-Hermite type. Assuming 
our Conjecture IB.ll we can apply the higher order variational equations for testing the 
integrability in the considered case. Calculations show that a logarithmic term appears in 
the second order variational equations, so the system is not integrable. 

If condition (|7.3p is satisfied, then either we obtain a system equivalent to (|7,4p . or to 
the following one 

q\ = -Xq 1 q 2 , q 2 = -ql A G C*. (7.10) 

If A belongs to one of iteme in table (|1.13|) (with k = 3), then the necessary conditions 
of Theorem 11.31 are satisfied. It seems that for all these values the force admits two or 
three independent polynomial first integrals. A family of forces with the above properties 
can be find for an arbitrary k > 2. We discuss these families in the last subsection of this 
section. 

At the beginning, we assumed that the force admits at least one Darboux point. Hence 
now, we consider all cases when there is no a Darboux point. By Theorem 16.31 we can 
assume that the force has the form given by (|6.37p with k = 3. If in (16.370 p = 2, then 
n\ = 1 or ri2 = 1 and the system is not integrable. When p = 1, m = 3 the force has the 
following components 

= (<?2 - z 1 qi)(q 2 - z 2 qi), F 2 = (q 2 - z 1 q 1 )(q 2 (2z 1 - z 2 ) - zfqi). (7.11) 

Let Qi = qi and Q 2 = q 2 — z\q\. Then in new variables the Newton equations with the 
above force read 

Qi = Q2{Q 2 -aQ 1 ), Q 2 = ~aQl a E C, (7.12) 
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and have one first integral 

I = \pI + \aQl P 2 = Q 2 . 

If a = 0, then system (|7,12p is a sup erintegr able Hamiltonian system with the Hamilton 
function ^ 

H = P1P2 - -<2 2 , 

and two additional first integrals 

h = P2, h = 4P 2 (QiP 2 - Q2P1) + Qi 

If a 7^ 0, we can assume that a = 1. For this case, using a direct method, we checked that 
the system (|7,12p does not admit a second polynomial first integral of degree lower than 
12 with respect to the momenta. 

7.2 Case k = 4 

To shorten the exposition, we omit proofs and many details of calculations in this sub- 
section. At first, we recall from paper that for k = 4 set of Aj = Aj — 1 with 
Aj satisfying assumptions of Theorem 11.31 has elements listed in Table 1, and X±i = for 
2 < I < 4. 



element 


{Ai,A 2 ,A 3 ,A 4 } 


1 


{-1,-1,2,2} 


2 


{-5/8,5,5,5} 


3 


{-5/8,2,20,20} 


4 


{-5/8,27/8,27/8, 135} 


5 


{-5/8,2,14,35} 



Table 1: Elements of set I4 4. 



7.2.1 Four Darboux points 

If we assume that a homogeneous force F = {Fx^F^) has four simple Darboux points, 
then, by a proper choice of coordinates, three of them can be located arbitrarily on the 
Riemann sphere. Thus, we can assume that one of them is at infinity, and two of them at 
z = and z = 1, respectively. Under these assumptions we have 



Fx = a n q{ + ai2<7i<?2 + ai3<7i<7 2 > ^ ^ 

F 2 = a 2 xq\q 2 + a 22 qxq 2 + a 23 q 2 , 



where 

an + ai 2 + ai3 = 0-21 + a 22 + a 23 = 1, 
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and 

011^23 / 0, and a u / a 2 i, and a i3 / a 2 3, and an - a 2 i / a i3 - a 23 . 
Now z-coordinates of Darboux points are 

z = 0, z± = 1, z 2 = , z 3 = CO. 

Ol3 — 0,23 

The respective non-trivial Kovalevskaya exponents at these points are the following 

a o\ . . (ai - a 3 )aia 3 a 3 
Ao = — , Ai = a 3 -ai, A 2 = 5- ■ 3, A 3 = -- 



an ana§ + ai 2 aia 3 + a i3 aj ' ai 3 + a 3 

where 

ai = a 2i - an, a 2 = a 22 - a i2 , a 3 = a 23 - a i3 . 

Now we can take the admissible values of Aj, i = 1, 2, 3, 4 from Table 1 and reconstruct 
Fi and F 2 . 

For (A , Ai, A 2 , A 3 ) = (-1,-1,2,2) we obtain 

F 1 = j[-(2a 2 + l)q\ + 5a 2 gig 2 + 3(2 - a 2 )gf], 

,? ( ^ ^ (7-14) 



F 2 = ^[(a 2 + 3) gi + (2 - 03)95], a 2 |-^,2 
for (Ao,Ai,A 2 ,A3) = (5,5,-5/8,5): 



i 7 ! = y[(«2 - 7)g 2 + 5a 2(M2 + 6(2 - a 2 )q 2 2 ], 

92, 



F 2 = y [6(a 2 - 7)qj + 5(9 - a 2 )gig 2 + (2 - a 2 )qj], a 2 {2, 7}, 
for (A ,A!,A 2 ,A 3 ) = (2,20,-5/8,20): 

F 1 = g-[10(a 2 - 22)g 2 + lla 2 g ig2 + 21(11 - a 2 )g 2 ], 

F 2 = ||[30(a 2 - 22)q 2 + ( 6 60 - 29a 2 )gig 2 + (11 - a 2 )q% a 2 {11,22}, 
for (A , A l5 A 2 , A 3 ) = (-5/8, 27/8, 27/8, 135): 

Fi = [27(22 - ha 2 )ql + 220a 2 q iq2 - 17(22 + 5a 2 )g 2 ], 

a? o 1 ( 22 221 

F 2 = Y^7jt 81 ( 22 " 5fl 2)9? + 4Wa 2 qiq2 - (22 + 5a 2 )g 2 ], a 2 — j 

and for (A , Ai, A 2 , A 3 ) = (2,35,-5/8, 14): 

F 1 = f^[7(2a 2 - 77)g 2 + 16a 2 q iq2 + 15(37 - 2a 2 )g|], 

F 2 = ||[21(2a 2 - 77)g 2 + 4(399 - 10a 2 )g ig2 + (37 - 2a 2 )gl], a 2 I —, — 
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(7.15) 



(7.16) 



(7.17) 



(7.18) 



Permutations of allowable values of (Ao, Ai, A2, A3) give equivalent forces in all the cases 
listed above. 

As it has already been mentioned, there is no forces satisfying conditions A.l and A. 2 
formulated before Theorem 16.11 with a smaller number of simple Darboux points. 

For forces (|7.14p - (l7.18p the necessary integrability conditions formulated in Theo- 
rem 11.31 are satisfied. Thus, to check if they are integrable, we applied the method of 
higher order variational equations. Our calculations for forces (|7.15p - (|7.18p show that 
the only integrable cases correspond to known Hamiltonian systems with homogeneous 
potentials. For a Newton equation with force (|7.14|) the variational equations do not re- 
duce to the product of two Lame equations and the application of higher order variational 
equations is an open problem. 

We also checked whether the Newton equations with forces (|7.14p ^ (l7.18p admit a single 
polynomial first integral for specific values of a 2 . To this end, we applied a direct method, 
and we looked for first integrals which have degree lower than 8 with respect to the 
momenta. The results of these calculations are following. 

Force (I7.14p admits a polynomial first integral only for a 2 = —3, but in this case the 
Newton equations are Hamiltonian and integrable. 

Force (|7.15|) has always a polynomial first integral 

h = Qi(li ~ Q2)q2 [(«2 - 7)qi - (a 2 - 2)q 2 ] 2 - 5(o 2 - 7)q 2 p\ - 5(a 2 - 2)q 1 pj 
+ 5[(a 2 - 7)gi + (a 2 - 2)q 2 ]p 1 p 2 . 

Furthermore, for a 2 G {—3,9/2, 12} the system is Hamiltonian and integrable. 

For (|7.16|) if a 2 = the system is Hamiltonian and integrable while for a 2 = i|(21±\/7) 
is Hamiltonian but non-integrable. Except these Hamiltonian cases only for a 2 = 22 
force (|7.16|) admits a polynomial first integral of the form 



Force (|7.17p admits polynomial first integrals only for a 2 = {0, ±||y^||^} when the 
system is Hamiltonian and non-integrable. 

Finally, force (|7.17|) admits a polynomial first integral only for a 2 G {0, ^| (285=f4\ // 86)} 
when the system is Hamiltonian and non-integrable. 

7.2.2 Smaller number of simple Darboux points 

If the considered force admits fewer than four simple Darboux points, than we have to 
investigate all the cases when condition A. 2 formulated before Theorem 16. II is not satisfied. 
It is easy to check that the case when k = 4, and the system possesses three Darboux 
points is impossible. 

Hence, we consider forces with two simple Darboux points. It can be shown that in 
this case the force is given by 



/ = I2pl + 8(71 {q 7 2 + Ap 2 2 qi) - 4 P2 (2 Pl + 3p 2 )qj + 3 9 | - Sqfql 




F\ = Qi(qi ~ Q2) [(A 2 - l)gi + (Ai - l)A 2 g 2 ], 
F 2 = q 2 (q± - q 2 )[(Xi - l)q 2 + (A 2 - l)Aigi]. 



(7.19) 
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provided that Ai 7^ 1 and A2 7^ 1. If Ai = 1 or A2 = 1, then the above force has infinitely 
many Darboux points and, as we showed in Section T6.41 the corresponding Newton system 
is integrable. The above force satisfies the necessary conditions of Theorem 11.31 iff Ai and 
A2 take values from table f)l . 13|) . Let us notice that to derive (|7.19p we have to use a 
generalized relation (16.230 as described in Remark 16.31 

Components of force (17.190 have a common factor q\ — q 2 , so we can use Theorem 16. 21 
Applying Proposition 16.21 we easily transform the force to the form Fi = q 2 Fi, where Fi 
are homogeneous polynomials of degree 2 for i = 1,2. It appears that deg^ F 2 = 2 iff 
Ai 7^ A2. Hence, by Theorem 16. 2[ the Newton system with force (|7.19p is not integrable if 
Ai + A 2 . 

If A = Ai = A2 7^ 1 is given by an item from table (|1.13|) . then, in some cases, we were 
able to find one first integral. For example, if A = 6, then the first integral has the form 

I = Q2?i + q\p\ - (qi + Q2)piP2 + 5qxq 2 (q 2 - qi) 3 - 

Also when Ai 7^ A2 we can find one first integral. For example, for Ai = 3/8 and A2 = 6 
the first integral reads 

I = P2(qm - Q2Pi) - ^qiq^qi - q2) 2 - 

Now we consider the case with only one Darboux point. It can be shown that in this 
case the force is given by 

Fx = aq± q 2 (qi + 92), F 2 = q 2 (qx + <?2)[(1 - a)qx + 52], ol 7^ 1, or (7.20) 
Fx = aq\q 2 , F 2 = q 2 (q 2 + aqx), or (7.21) 

Fx = Xqxql F 2 = q*. (7.22) 

By Theorem 16. 2} the Newton system with force (|7.20p is not integrable. For a = 
the Newton equations with force (|7.2ip are Hamiltonian and integrablejwhile for a 7^ 
they are also Hamiltonian, however nonintegrable, as it was proved in [22|. The case of 
force (|7.22p is discussed in the last subsection of this section. 

Assuming that Conjecture IB . 1 1 is valid, and applying the higher order variational equa- 



tions method similarly as it was done for Hamiltonian systems in 20|;|22|, it is possible to 
prove the following. 

Lemma 7.1. The only integrable Newton systems with k = 4 possessing a multiple Dar- 
boux point are equivalent to the following one 

q\ = -qx{axxq\ + qi), h, = -q^iaxxqf + qi), 

with first integrals 

h = q\P2 - q2.P1, h = Tjiaiipj +pl) + -^(axxql + qi) 2 - 
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Finally, let us consider the integrability problem for forces without a Darboux point. 
Then by Theorem 16,31 the only forces that can admit two first integrals are either of the 
following: 



Fx 


= 92 [92 (91 + aq 2 ) +q\], 


F 2 


= 92(9i +92), 


or 


(7.23) 


F x 


= o, 


F 2 


= -9192) 


or 


(7.24) 


Fx 


= 9192, 


F 2 


= -"9192 , 


or 


(7.25) 


Fx 


= o, 


F 2 


- a 3 

— 9i> 


or 


(7.26) 


Fx 


= 9l(9l + 92), 


F 2 


- a 3 

— q 2 - 




(7.27) 



The force field (|7.23|> for any a possesses a first integral that is a nondegenerate quadratic 
form with respect to the momenta with constant coefficients, thus it is always Hamilto- 
nian, however we found the second functionally independent first integral only for a = 0. 
Calculations were made up to degree 8 in the momenta. The force with components (|7,24p 
has one first integral I = px and there is no the second polynomial first integral up to 
degree 9 in the momenta. Newton system governed by (I7.26P is integrable and Hamilto- 
nian. The last force (|7.27|) is non-Hamiltonian with one first integral / = 2p| + q\ and the 
direct method does not yield the second polynomial first integral up to degree 10 in the 
momenta. 

The most puzzling are forces (IT.25I) . The application of the direct method up to 
degree 7 in the momenta yields one integrable Hamiltonian case for a = — 1, the another 
one integrable non-Hamiltonian for a = —11/6 with two first integrals which are both of 
degree four in the momenta 

h = 2A PlP 3 2 + 3q1(Aqlp 2 2 + l2 qi q 2 pxP2 ~ MIpI) + 16g?gi, 

I 2 = 162p 3 1 (q 1 p 2 - q 2 px) ~ 9^(4^^ - 20qxq 2 pxp 2 + IZqjpj) + 16qfql 

Moreover we found that for some rational and irrational values of a one polynomial first in- 
tegral exists. This is the case for a G {0, -7/2, -4, -5, -10, -11, (-38+21>/3)/ll, -187/7}. 
For these values of a first integrals are quite complicated as it is shown by some examples 
below: 

I = 2px(qxp 2 - q 2 px) + qfq 2 , 
I = 4pl-ql(q lP 2-4q2Px), 

I = Qp\p 2 + q\(q\p\ - 4qxq 2 pxp 2 + 15^p?) + 2q\ql, 

I = 176(38 + 2lV3)p?p| - 33 [( 39 + 2 3V3)9ip| - 4(39 + 23v / 3)9i9 2 PiP2 

- 10(7 + ^)q\q\p\p\ + 4(-3 + h^)qxq\p\p 2 + (3 - 5^3)?^] 
+ 264^ 9 |[3(2 + V3)q 2 lP 2 2 - 2q x q 2 pxp 2 - 3(-2 + V$)qlp\\ 

- 32(-7 + sV3)qfql 



a 

a 
a 

a 



7 

~2' 
-4, 

-5, 

-38 + 21^/3 
11 ' 
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7.3 Case k > 4 



We performed an analysis of the case with k = 5 in the same way as for the case with 
k = 4 in the previous section. It does not give any intersting examples of Newton's 
equations except for a case when the force admits only one simple Darboux point (see the 
last subsection) or when they do not admit any Darboux point. It should be mentioned 
that the bigger number of parameters make the analysis more complicated. 

Below we describe a few general facts which are valid for arbitrary k > 5. Let us 
consider a case when the force has exactly k Darboux points. Then for each k set & 
contains two elements 



^±- ^fc + l,...,fc + J } . (7.2S) 

fe— 1 times 



fc— 2 times 

Let us note that for k = 14, 17, 19, 26, . . . this set is bigger, see [iij. We show the following 
fact. 

Lemma 7.2. Let F = (F\, F2) be a homogeneous force possesing k Darboux points, where 

(2) 

k = 1 + deg Fi. Assume that (Ai, . . . , Afc) = T) I, where Aj is i/ie non-trivial Kovalevskaya 
exponent of F at i-th Darboux point. Then F is equivalent to the following force 

k-l 

F 2 = ^(A ; + l+i)^- 1 - i ^. 

Proof. We can assume that one Darboux point is at infinity. The polynomial g(z) has 
degree k — 1 and the degree of polynomial is smaller than k — 1, see formulae (|6.5p . 
We put 

fc-i fc-i 



^(z) = 9i J = g k _ x ]J(z- Zi ). (7.30) 



i=0 i=l 



Then 



fe(z) _ ^ g^) 1 _ y 1 1 
g{z) h(zi) z- Zi ^AiZ-Zi 



Since Ai = • • • = A^-i = k + 1, we have 

jfc-l 



As result 



«=i x 7 
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Direct substitution of the explicit form of g{z) yields 



fc-i 

A = sxi5> + 1 + 



and homogenization of (f\, f 2 ) and multiplication by k + 1 finishes the proof. □ 
One can check directly that 

k 

I=(k+ l)pi{qxp 2 ~ q 2 Pl) +Y J ^-iQi +2 ~ i Q2~\ 

i=l 

is a first integral of Newton's equations with force (|7.29p . For this force the necessary 
integrability conditions of Theorem 11,31 Morales-Ramis theory are satisfied. Thus we can 
only apply a direct method for searching the second first integral. Here the problem is 
not trivial as force (|7.29p depends on many parameters. We have not been able to make a 
general search for first integrals of degree higher than 2 with respect to momenta. In all 
cases when a first integral was found, the system appeared to be Hamiltonian or trivially 
integrable. 

Lemma 7.3. Let F = (F\, F^) be a homogeneous force possesing k Darboux points, where 
k = 1 + deg -Fj. Assume that (Ai, . . . , A&) = XpjL where Aj is the non-trivial Kovalevskaya 
exponent of F at i-th Darboux point. Then F is equivalent to the following force 

k-l k-X 

Fl = J2(i + l- fc)<?i<tf-V 2 -\ F 2 = £(i - Vmt^qi (7.31) 
i=l i=l 

Proof. We can assume that one Darboux point is at the infinity and for it Ao = —1. The 
other with Afc_i = A_ = —1 can be locate at z^-i = 0. For Darboux points Z\, . . . , -Zfc_2 
we have Ai = • • • = A^-2 = A = k — 2. Thus 

h _ 1 v 1 1 1 - 1 V 1 I 1 1 

g A ' y. — ■ \ - — ■ i • \ ^-^ - — - 1 \ A \ 
or, in other words 



5 A^z-Zj A_z-z fc _i A^z-Zj \A_ Ajz-Zk-i 

1=1 1=1 



A \A_ A/ z-zn fc-2 (k — 2)(z — Zk-i) 

We recall that fi = h and $2 = 9 + zh. Thus, putting z^_\ = and multiplying by k — 2, 
we obtain 

Substitution of ()7.30p into the above expressions and homogenization gives (|7.3ip . □ 



41 



For force (I7,31|) the necessary integrability conditions of Theorem 11.31 are satisfied. 
Hence, as in the previous case, we can only apply a direct search for two first integrals, 
however conclusions from such an investigation are the same as for force (|7.29p . 



7.4 An integrable family of Newton equations 

Here we consider a force given by 

F x = Xqiqt 2 , Fl = A G C*, k > 2. (7.32) 

This force has one Darboux point at infinity, and it satisfies the necessary conditions 
for integrability of Theorem 11.31 provided A belongs to an item of table 11.131 It appears 
that for all these values of A, the Newton equations are really integrable. Moreover, this 
integrability is highly non-trivial. 

Force (|7.32p admits one obvious first integral 

h = \pl + \ql (7.33) 

Assume that 

k 

X = p + -p(p-l), p€Z*. (7.34) 

Then for each k and p there exists an additional first integral of degree p with respect to 
the momenta for p > 0, for p < this degree is 1 — p. For example, if k = 3 and p = 6 
this first integral reads 

h = 28pl(q 1 p 2 - q 2 Pi) - 70gfp|(6gip 2 ~ I2P1) + I5q%p 2 (21qip 2 - qzpi) - 15gig|, 
while for k = 6 and p = —5 it takes the form 

h = llpip% + ll<?|p!(17giP2 - 2q 2 pi) - 3g| 1 p 2 (34gip 2 - + 3qiq\ 7 ■ 

We checked the above claim finding explicit forms of integrals for k < 10 and \p\ < 10. Let 
us note that for natural Hamiltonian systems with a homogeneous potential in all known 
integrable cases the additional first integral is of the degree smaller than 5 with respect to 
the momenta. 

If we assume that 



— + kp(p + 1) 



pGNU{0}, (7.35) 



the force admits two additional first integrals, so the corresponding Newton equations are 
superintegrable with three polynomial first integrals. Degrees of two additional integrals 
with respect to the momenta depend on p and are unbounded with respect to p. For 
example, for k = 3 and p = the additional first integrals read 

h = 6pi(<7iP2 - Q2Pl) + qlqh h = 6p? - q\{q\P2 ~ 3<? 2 pi), 
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while for k = 4 and p = 1 they have the following forms 

h = 240p lP l(q lP2 - q 2Pl ) + 2q\{\lbq\p\ - 2Q0 qi q 2Pl p 2 + ±qbl) ~ 22hq\q\, 
h = 10368^ + 16(625^ - 2500qlq 2 p lP 3 2 + 31^q\q\p\p\ + 1280 qi qip 3 lP2 + 4q$p$) 
- 1350(^(25^ - 80 qi q 2 pip 2 - 8qlp\) + 50625^- 

We checked the above claim as the previous one, finding explicit forms of first integrals. 
In all cases we checked that three first integrals I\, I 2 and /3 are functionally independent. 

For k > 5 the considered above two choices for A are the only admitted by Theorem ll.31 
For k < 5 we have additional admissible choices for A. For example, for k = 3 the four 
additional families depending on integer parameter p appear. Below we show first integrals 
for a few elements of these families: for A = — + |(1 + 3p) 2 

p = 0, I 2 = 192pi - lQqlpi{q\p 2 - 3q 2 pi) - q\q\, 

h = l92pl(q lP2 - q 2Pl ) + 2qj{qjp 2 - A qi q 2 p lP2 + 2Aq 2 p 2 ) + qfq 3 , 
p = 1, I 2 = 5184pfp| - 16pi(343^p| - W29q 2 q 2 p lP 2 2 - 672q iq 2 p 2 p 2 + %q\p\) 

- U7q 2 q 2 2 {49qjp 2 2 - 224 qi q 2Pl p 2 - 32q 2 p\) + 14406^, 

I 3 = 36288plp 4 2 {q lP2 - q 2Pl ) + 14p|(343gfpl - W2q 3 q 2Pl p 3 + 12264g 2 q 2 p 2 p 2 

- 12928gig|pfp 2 + 640g|pf ) - g|(7203gfp^ - 208544g? q 2 p\p\ + 272832g 2 gfp 2 p 2 . 

- 16128gigfp?p 2 - 256g|pf) + 588g 2 ^(147g 2 ^ - 280qiq 2Pl p 2 + IGq^pf) - 28812^, 

for A = -^ + ^(l + 4p) 2 

p = 0, / 2 = 294912p? + 15360g 2 p?(-gip 2 + 3q 2Pl ) - 16qf(q 2 p 2 2 - 6qiq 2P i P2 + 90g|p 2 ) - 9qfq 2 , 
I 3 = 28311552p[(g 1 p 2 - q 2 p x ) + bim9Qq 2 p\{q 2 p 2 - iqiq 2Pl p 2 + 16g 2 p 2 
+ 1536qi Pl (-qfpl + lq\q 2 pxp\ + 14giglp 2 p2 + 84g 2 i p?) - 8g?gl(5g 2 ^ + \32q x q 2 p x p 2 

- \W8q\p\) - 21q\q\. 

For A = — 2j + ^(1 + 5p) 2 and for A = — ^ + ^j(2 + bp) 2 expressions for I 2 and 1$ are 
more complicated and we do not write them here. 

Our calculations show that for additional families of A's there always exist two addi- 
tional functionally independent first integrals. 

The investigations shortly presented above allowed us to formulate Conjecture 11.11 
Without doubt, the integrability analysis of equations (|1 . 17f) presented above is not satis- 
factory, and a separate paper [23j will be devoted to more involved investigations. 

Appendix A 

Here we characterize in two ways forces that are potential in the sense of Definition [2j 
Lemma A. 4. Assume that Newton's equations 

q = p, P=-F(q), q,peC n , (A.l) 
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admit a first integral of the form 



I=-p T Kp+W(q), 



(A.2) 



where K G GL(n, C), K = K T . Then force F is potential and equations (lA.lj) are 
Hamiltonian with respect to a symplectic structure given by the following matrix 



J 



K~ l 
-K' 1 



(A.3) 



and I becomes the Hamiltonian function. 

Proof. Since I is a first integral of equations (jA.lh . we have 



V 



dq 



+ p KF = 0, 



so 



oq oq 



Hence 



p 




-F 





K- 1 
-K- 1 



dl 
dq 

m 

dp 



In this way we proved the second claim of our lemma. To prove the first one, let us 
notice that for every non-singular matrix K there exists a matrix A £ GL(n, C) such that 
K = AA T . Hence, in new variables 



Q = Aq, P = Ap, where AA 1 



equations (|A.1|) read 

Q = P, P = -W(Q), V(Q) := W(A~ 1 Q), 
so they are canonical Hamilton's equations with 

H := I{A~ l Q, A- l P) = X -P T P + V(Q), 

as the Hamiltonian function. 



(A.4) 
(A.5) 

(A.6) 
□ 



The second characterization of potential forces is the following. If a force F(q) is 
potential, then, by Definition [21 there exists a non-singular matrix A such that 



F A (Q) = A~ 1 F(AQ) 



dV(Q) 
dQ ' 
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for some scalar function V. Hence, putting q = AQ we obtain 

KF (g) = ^M, K~ 1 = A T A. (A.7) 

This shows that F is a potential iff there exists a non-singular matrix K such that differ- 
ential form 

n n 

w = (KF, dg) :=J2J2 K H F J d 9i, ( A - 8 ) 
i=\ j=i 

is exact. 

The condition dw = is equivalent to the following set of equations 

i=i 



K pi —^ - K ri —± = 0, 1 < p < r < n. (A.9) 
y dq r dq P/ ' 



For a given force in the above equations the elements of matrix are the unknowns K = [Kij] 
with i < j. Since we consider the polynomial force and degFj = k — 1, all n(n — l)/2 
equations (|A.9h are polynomial and homogeneous of degree k — 2, Hence, they give rise to 
a system of homogeneous linear equations for unknowns . 
As an example we consider the case n = 2 and k = 3. We put 

T? °11 2 , , a 13 2 Tp a 21 2 , , a 23 2 /. ,m 

F\ = —q 1 +ai 2 qm + ~Y<1%, F 2 = —q x + a 22 qiq2 + ^q 2 - (A.10) 
Then equations (jA.9h are equivalent to the following linear system 

Lk = 0, L= ( ai3 ° 2 3-«i2 -«2 2 A (An) 

V«12 «22-ail -021/ 

and k = [Kn, K\ 2 , K 22 ] T . This system has always a non-zero solution. If one of 2 x 2 
minors Lj of matrix L 

Li = ai 2 (ai 2 -a 23 )+ai 3 (a 22 -an), L 2 = a 2 i(ai 2 -a 23 )+a 22 (a 22 -an) , L 3 = ai 2 a 22 -ai 3 a 2 i, 
is non-zero, then it has the follwing solution 

Kn = a 2 i(ai 2 - a 23 ) + 022(022 - an), K i2 = a n a 2 i - ai 2 a 22 , 
K 22 = ai 3 (a 22 - an) + au(a 12 - a 23 ), 
and matrix K given by this solution is not singular iff 

Oi2[«i2a2i - aii022 - 2a 2 ia 2 3] + 012 [(an - 022)022023 - 013021(011 - 3a 2 2) + a 2 ia| 3 ] 
+ Oi 3 [(an - 022) ((on - 022)022 + 021023) - 0130^] 7^ 0. 

(A.12) 

For values of ajj for which all minors Lj vanish one can find a two-parameter family 
of solutions of (jA.lip and always obtain a non-singular matrix K. 

For n = 2 and k > 3, as well as, for n > 2 equations (|A.9h give rise to an overdetermined 
system of linear homogeneous equations, so a non-zero solution exists only if the rank of 
matrix L is not maximal. As entries of L depend on coefficients of polynomials F\ and 
F2, this gives very restrictive conditions for a force F = (F\,F 2 ) to be potential. 
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Appendix B 



Assume that we know a particular solution <p(t) of a holomorphic system 

—x = v(x), xeC, tec. (B.l) 

Then appart from variational equations 

^ = Am, A{t) = j£(fp(t)), (B.2) 
one can use higher order variational equations. The idea is following. Let us put 

where e is a formal small parameter. Inserting the above expansion into equation (|B. 1|) 
and comparing terms of the same order with respect to e, we obtain the following chain 
of linear non-homogeneous equations 

l$< fc > = A(t)£W + f k (^\. . . , Z^-V), k = 1, 2, ... , (B.3) 

where / 1 = 0. For a given k equation (jB.3[) is called the A;-th order variational equation. 
We denote by X(t) the fundamental matrix of the homogeneous system, i.e., nxn matrix 
satisfying 

^X = A(t)X, X(0) = E, 

where E is the identity matrix. Then solutions of k-th order variational equations for 
k > 1 are given by 

= X(t)c(t), (B.4) 

where c(t) is a solution of 

±-c = X-\t)f k . (B.5) 

There exists an appropriate framework allowing to define the differential Galois group 
of the k-ih order variational equation, for details see 24j; |27|. For Hamiltonian systems 
the following theorem was proved in 27]. 

Theorem B.l. Assume that a Hamiltonian system is meromorphically integrable in the 
Liouville sense in a neighborhood of a phase curve V which is not an equilibrium point. 
Then the identity component of the differential Galois group of the k-th order variational 
equations is Abelian for any KN. 

The main problem with an application of the above theorem lies in the fact that it 
is very difficult to investigate the differential Galois group of higher order variational 
equations. However, if the first order variational equations are a direct product of Lame 
equations, then there is a local criterion which allows to decide whethere this group is 
Abelian, for details and practical applications see 24; 27; 2C; Izj]. 
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Let us note that the main construction presented in [27f | is valid also for non-Hamiltonian 
systems. Thus, we can also use it for the Newton's equations. Unfortunately, there are 
many theoretical problems in extending Theorem IB.ll in such way that it should apply 
to the Newton equations. Nevertheless, there are many facts strongly suggesting that the 
following conjecture is true. 

Conjecture B.l. Assume that a Newton system of the form (II. 5p is meromorphically 
integrable in the Jacobi sense in a neighborhood of a phase curve V corresponding to a 
Darboux point. Then the identity component of the differential Galois group of the k-th 
order variational equations is Abelian for any k E N. 
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